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ratio  of  specific  heats,  1.4  for  air 

boundary  layer  thickness,  defined  by  location  where  M=0.95Me 
reference  boundary  layer  thickness  at  M=2.8  (6.1mm) 
reference  boundary  layer  thickness  at  M=4.9  (10.2mm) 
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CHARACTERIZATION  OF  THE 
FLOWFIELD  NEAR  A  WRAP-AROUND  FIN 
AT  SUPERSONIC  SPEEDS 

1.  Introduction 

1.1  Wrap-Around  Fins 

Wrap-around  fins  (WAFs)  have  been  used  by  designers  for  several  years  on  low- 
speed  tube  launched  missiles  and  dispenser-launched  sub-munitions.  The  term  “wrap¬ 
around  fin”  usually  refers  to  a  projectile  stabilizing  or  control  surface,  which  has  the  same 
curvature  as  the  missile  body,  and  can  be  wrapped  around  the  projectile  until  deployment 
(see  Figure  1.1).  Since  stealth  capability  has  become  a  design  parameter  for  many  aircraft, 
WAFs  have  become  even  more  attractive  for  their  reduced  cross-section  and  stowability. 
Wrap-around  fins  can  also  simplify  the  design  of  airframes  that  integrate  the  weapon  in 
partial  concealment,  avoiding  complications  associated  with  fin-fuselage  contact. 

While  WAFs  enable  several  design  possibilities,  several  stability  anomalies  are  in¬ 
herent  for  missiles  employing  them!1'  6>  52’ 113’  1191,  the  most  recognized  of  which  is  a  roll 
reversal  observed  near  transonic  conditions.  Also  due  to  the  asymmetric  fin  geometry, 
missiles  with  WAFs  display  a  pitch-yaw  coupling  at  all  speeds  not  present  on  conventional 
missiles.  During  recent  ballistic  range  tests!1’  1131,  a  second  rolling  moment  reversal  was  ob¬ 
served  over  a  small  range  of  high  supersonic  speeds  (M«4.5-4.7)  on  a  WAF  configuration. 
In  this  regime,  yawing  moment  reversals  were  also  detected.  Vitale  and  Abate I11^  have 
proposed  that  this  loss  of  static  stability  may  be  related  to  the  complex  shock  structure  in 
the  fin  region.  Cross-flow  induced  by  missile  pitch  and  spinning  may  also  be  a  contributing 
factor.  Interaction  with  the  missile  bow  shock  is  also  plausible  at  high  pitch  angles. 

The  majority  of  WAF  experiments  have  focused  on  ascertaining  stability  character¬ 
istics  via  sub-scale  flight  tests!4'  5. 6,  52, 104, 113,  H7]?  most  of  which  emphasized  the  subsonic 
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and  transonic  flight  regimes!4,  s,  6,  52,  H7],  These  experiments  provide  no  detailed  flowfield 
measurements,  and  do  little  to  promote  an  understanding  of  the  flowfield. 

A  small  number  of  numerical  simulations  have  been  performed  on  missiles  employing 
wrap-around  fins.  Some  have  focused  on  the  design  of  lift  and  control  surfaces  by  simulating 
the  fin  alonet58h  However,  the  majority  of  these  simulations  have  focused  on  characterizing 
the  structure  of  the  fin  shocks  and  their  interaction.  The  Euler  equations  were  solved  by 
Vitale,  et  all113!  on  a  swept  WAF  configuration  with  no  fin  thickness.  The  Euler  equations 
have  also  been  solved  by  Abate  and  Cook!3!  for  subsonic,  transonic,  and  supersonic  flow 
conditions  on  a  WAF  configuration  using  unswept  10%  bi-convex  airfoils.  While  such 
inviscid  simulations  have  been  able  to  predict  the  existence  of  the  sonic  rolling  moment 
reversal,  they  have  failed  to  predict  the  second  rolling  moment  reversal  at  high  Mach 
numbers.  More  recently,  Edge!28!  has  solved  the  laminar  Navier-Stokes  equations  over  a 
WAF  geometry  with  a  relatively  low  fin  aspect  ratio  and  a  round  leading  edge  for  Mach 
numbers  ranging  from  1.3  to  3.0.  To  date,  no  WAF  geometry  has  been  simulated  using 
turbulence  models. 

1.2  Objectives  of  Present  Research 

The  primary  goal  of  the  present  research  was  to  characterize  the  flow  structure  near 
WAFs  at  high  supersonic  speeds.  This  would  be  accomplished  through  a  systematic  nu¬ 
merical  and  experimental  analysis  of  the  mean  and  turbulent  flowfield  in  the  vicinity  of  a 
single  WAF.  An  understanding  of  the  flowfield  near  WAFs  is  critical  to  further  development 
of  such  configurations,  given  the  dependence  of  stability  characteristics  on  Mach  number. 
Determination  of  the  flow  structure  near  a  single  non-spinning  WAF  is  an  essential  first 
step  toward  this  understanding. 

A  second  objective  was  to  contribute  a  complete  set  of  mean  flow  and  turbulence 
data  on  this  shock/boundary-layer  interaction  flowfield  for  numerical  turbulence  model 
validation.  The  data  gathering  criteria  of  Settles  and  Dodson!96!  has  been  used  for  guidance. 
Prior  to  this  study,  no  detailed  flowfield  measurements  (mean  flow  or  turbulence)  existed 
for  WAF  missile  configurations. 
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Finally,  the  suitability  of  a  widely  used  commercial  simulation  package,  the  General 
Aerodynamic  Simulation  Program  (GASP) I7!,  with  the  turbulence  model  of  Baldwin  and 
Lomax!11!  was  to  be  investigated  for  this  flowfield. 

1.3  Overview  of  Present  Research 

As  a  first  step  toward  understanding  the  flow  structure  near  WAFs,  a  simple  model 
has  been  investigated  which  consists  of  a  single  wrap-around  fin  mounted  on  a  partial  fuse¬ 
lage  (Figure  3.1).  This  simplified  model  allowed  experimental  data  to  be  obtained  at  much 
higher  resolution  than  would  have  been  possible  on  a  full-body  four-finned  configuration 
scaled  to  fit  in  the  available  tunnel  space.  Also,  the  single- WAF  static  model  isolates  the 
effects  of  fin  curvature  from  the  effects  of  upstream  cross-flow  and  interaction  of  the  multi¬ 
ple  fin  shocks.  The  flow  around  this  single- WAF  configuration  has  been  investigated  both 
numerically  and  experimentally,  with  the  objective  of  quantifying  the  flow  structure  in  the 
region  near  the  fin/body  juncture.  Numerical  and  experimental  results  are  examined  and 
compared  to  gain  an  understanding  of  the  flow-field  in  the  vicinity  of  a  WAF. 

Measured  flowfield  data  were  obtained  in  the  AFIT  Mach  3  wind  tunnel  and  in  the 
new  AFIT  Mach  5  wind  tunnel  on  a  model  composed  of  a  single  wrap-around  fin  on  a 
ceiling-mounted  semi- cylindrical  fuselage.  Shadowgraph  and  schlieren  photographs  were 
also  obtained  for  flow  visualization.  In  the  Mach  3  wind  tunnel,  the  flow  around  the 
test  article  was  surveyed  at  several  stations  along  its  length,  concentrating  on  the  region 
near  the  fin.  The  flowfield  was  also  explored  with  the  model  in  wall-mounted  configura¬ 
tions  to  enable  measurements  closer  to  the  fin  surface.  Detailed  mean  flow  measurements 
were  obtained  using  conventional  cone-static  and  pitot  pressure  probes,  as  well  as  hot- 
film  cross-wire  probes.  The  flow  data  measured  with  the  hot-film  probes  were  also  used 
to  estimate  turbulence  quantities  using  the  hot-film  anemometry  methods  outlined  by 
Bowersox  and  Schetzt14’  16>  16)_  These  techniques  have  been  used  recently  in  several  ex¬ 
periments!26’  74'  75’  94J  to  quantify  the  effects  of  pressure  gradients  on  turbulent  boundary 
layers.  In  the  Mach  5  tunnel,  mean  flow  (pressure)  measurements,  as  well  as  shadowgraph 
and  schlieren  photography  were  obtained. 


1-4 


Experimental  data  were  compared  to  numerical  results  obtained  with  a  widely  used 
numerical  simulation  package  (GASP  V3.0617!)  employing  the  Baldwin-Lomaxt11!  algebraic 
turbulence  model.  Taken  in  concert,  the  experimental  and  numerical  information  are 
examined  with  a  view  toward  characterizing  the  net  effect  of  the  complex  flowfield  in  the 
vicinity  of  the  WAF  on  aerodynamic  loading. 

1.4  Outline  of  Document 

This  section  provides  a  road  map  of  the  chapters  that  follow.  First,  Chapter  2  pro¬ 
vides  a  survey  of  research  that  has  been  directed  toward  the  solution  of  other  shock/boundary- 
layer  interaction  flows  sharing  some  characteristics  of  the  flowfield  near  supersonic  wrap¬ 
around  fins.  The  experimental  methods  used  in  this  investigation  are  outlined  in  Chapter  3. 
Details  relating  to  the  reduction  of  hot-film  anemometry  data  are  presented  in  Appendix  B. 
The  uncertainty  of  these  measurements  and  an  analysis  of  the  effect  of  these  uncertain¬ 
ties  on  the  presented  results  are  provided  in  Appendix  C.  Chapter  3  also  details  the 
facility,  apparatus,  instrumentation,  and  model  used  to  obtain  the  experimental  results 
at  Mach  3.  A  more  detailed  description  of  the  Mach  5  wind  tunnel  constructed  for  this 
research  is  included  in  Appendix  D.  The  numerical  strategy  used  is  outlined  in  Chapter  4, 
and  the  governing  equations  and  turbulence  model  that  were  employed  are  detailed  in 
Appendix  A.  The  results  of  the  experimental  and  numerical  investigations  at  Mach  2.8, 
including  an  analysis  of  the  flowfield  near  the  fin,  are  presented  in  Chapter  5,  while  the 
results  at  Mach  4.9  are  presented  in  Chapter  6.  A  set  of  conclusions  and  recommendations 
for  further  study  are  given  in  Chapter  7. 
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2.  Background 


In  the  mid-1960’s  research  on  shock  interactions  was  motivated  by  a  desire  to  quan¬ 
tify  the  extreme  heating  loads  imparted  on  bluff  bodies  by  single  and  multiple  impinging 
shock  waves.  Heirs  and  Loubsky  t39l  investigated  the  effects  of  shock  impingement  on  a 
cylindrical  leading  edge  in  a  shock  tunnel  at  supersonic  speeds,  though  the  experimental 
tools  of  the  day  limited  the  accuracy  of  their  heating  rate  measurements.  Meanwhile,  Ed- 
ney(29I  was  identifying  and  categorizing  the  flowfield  structures  produced  by  shock/shock 
interactions  for  several  shock-generating  geometries.  The  need  for  improved  high-speed 
bluff  body  aerodynamic  and  heating  prediction  methods  for  the  Shuttle  project  led  to 
several  experimental  investigations.  Holdent41-  8°1  provides  a  historical  review  of  several  of 
these  experiments  involving  single  and  multiple  shock  interactions  as  they  effect  aerother- 
mal  heating  loads.  This  chapter  surveys  the  experimental  and  numerical  investigations  on 
configurations  that  exhibit  similar  flowfield  characteristics  to  those  observed  in  the  vicinity 
of  supersonic  wrap-around  fins,  beginning  with  a  discussion  of  blunt  unswept  fins. 

2.1  Blunt  Unswept  Fins 

Over  the  past  three  decades,  several  efforts  have  been  directed  toward  understanding 
and  accurately  predicting  flows  having  shock/boundary-layer  interactions.  Efforts  have 
included  experimental  explorations  to  obtain  a  detailed  knowledge  of  the  flowfield,  as 
well  as  computational  studies  aimed  at  developing  and  validating  numerical  methods  and 
turbulence  models  for  such  flows. 

The  present  research  has  been  directed  toward  detailing  the  flow  structure  in  the 
vicinity  of  a  curved  fin.  It  was  anticipated  that  this  flow  structure  would  be  fundamen¬ 
tally  different  from  the  oft-studied  flow  around  blunt  unswept  fins,  given  the  significant 
cross-flow  component  of  the  mean  flow  induced  by  the  geometric  asymmetry.  However, 
the  flowfields  over  the  two  geometries  share  several  common  attributes,  and  an  under¬ 
standing  of  the  more  fundamental  flows  over  straight-fins  aids  interpretation  of  the  flow 
characteristics  near  wrap-around  fins. 
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Figure  2.1  Sketch  of  flow  interference  near  a  blunt  fin  (Hung  and  Kordulla,  1984). 

One  commonly  studied  geometry  of  particular  relevance  to  the  current  research  is  a 
blunt  unswept  fin  of  constant  thickness  and  a  rounded  leading  edge  mounted  on  a  flat  plate 
(Figure  2.1).  The  bow  shock  produced  by  the  blunt  fin  forces  the  incoming  plate  boundary 
layer  to  separate.  This  separation  creates  an  oblique  shock  which  interacts  with  the  bow 
shock  at  a  location  away  from  the  fin  root  (Figure  2.2).  A  supersonic  jet  forms  due  to 
the  Edneyl29)  Type  IV  interference  associated  with  the  interaction  of  the  separation  shock 
and  the  bow  shock.  The  entire  structure  is  referred  to  as  the  ‘A-shock’  pattern.  This 
shock  pattern  and  the  associated  vortical  flowfield  are  also  produced  ahead  of  a  cylinder 
extending  from  a  flat  plate  in  supersonic  flow. 

2.1.1  Experimental  Investigations.  As  early  as  1967,  Price  and  Stallings  !871 
investigated  supersonic  turbulent  separated  flows  in  the  vicinity  of  fin-type  protuberances. 
In  the  decade  following,  several  other  experimental  investigations  on  blunt  fins  in  supersonic 
flows  were  conducted,  most  notably  by  Kaufman,  et  al.t5^,  Sedney  and  Kitchens!92'93!, 
Dolling  and  BogdonofP24) ,  Ozcan!82'  83l,  Saida  and  Hattori!91],  Fomison!33'  34l,  and  Settles!95]. 

Sedney  and  Kitchens!92'93!  experimentally  studied  the  flow  near  cylinders  mounted 
on  a  flat  plate.  The  cylinder  height  and  diameter  were  varied,  as  well  as  Mach  number 
(1.5<M<4.5)  and  unit  Reynolds  number  (2xl06m-1<Re/£<19.3xl06m-1).  They  found 
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that  2,  4,  or  even  6  vortices  could  evolve  around  the  juncture  at  a  given  Mach  number, 
depending  on  the  unit  Reynolds  number.  The  number  of  vortex  pairs  generally  decreased 
with  increasing  Reynolds  number.  The  details  of  how  this  structure  varies  with  Re/£  is 
different  for  each  Mach  number.  However,  the  unit  Reynolds  number  (and  hence  number 
of  vortices)  was  found  to  have  little  effect  on  the  location  of  the  primary  separation  or 
attachment  lines  in  front  of  the  obstacle. 

Dolling  and  Bogdonoffl24!  demonstrated  (using  data  from  several  experiments!50’  87, 1201  j 
that  over  a  large  range  of  Mach  number  and  incoming  boundary  layer  thicknesses  {S/D), 
the  leading  edge  diameter  ( D )  was  a  suitable  scaling  parameter  for  the  centerline  plate 
pressure  ahead  of  the  fin.  The  upstream  influence  was  always  found  to  be  between  two 
and  three  diameters  upstream.  They  also  showed  that  for  a  given  Mach  number  (M=3, 
Re/£=65xl06m-1±5%),  the  entire  flat  plate  surface  pressure  is  scaled  by  the  leading  edge 
diameter.!24!  This  correlation  was  later  corroborated  by  Fomison.!33- 34!  It  was  found  that 
even  the  leading  edge  surface  pressure  ratios  were  correlated  with  y/D  for  a  given  Mach 
number,  provided  that  the  shock  wave  structure  was  well  clear  of  the  boundary  layer 
(i.e.,  S/D  <  4).  In  essence,  Dolling  and  Bogdonoff  showed  that  the  spatial  extent  of  the 
interaction  near  a  straight  unswept  fin  is  dominated  by  inviscid  characteristics  of  the  flow. 
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FomisonI33'  341,  investigating  blunt  unswept  fins  with  a  semi-cylindrical  leading  edge 
extending  from  a  flat  plate,  detailed  the  effect  of  fin  thickness  and  incidence  on  the  flow 
(Moo=2.4,  Re/£=2.6xl06m-1).  His  oil  flow  visualizations  showed  that  both  thick  and 
thin  fins  produced  the  same  number  of  vortices  (4  at  this  relatively  low  Re/^).  This  was 
consistent  with  the  hypothesis  of  Sedney  and  Kitchens!92'  93I  that  the  number  of  vortices  was 
driven  by  the  unit  Reynolds  number,  and  not  obstacle  size.  Fomison’s  oil  flow  experiments 
also  demonstrated  that  fin  incidence  can  have  a  dramatic  effect  on  the  flow  topology.  It  was 
conjectured  that  at  sufficiently  high  incidence  angles,  a  secondary  separation  would  appear 
on  the  wall.  This  has  been  shown  to  occur  for  sharp  swept  fins  at  high  incidence  angles.!33! 
Stolleryl80’  1031  has  reviewed  several  experiments  which  generate  flows  involving  glancing 
shock/boundary  layer  interactions,  including  blunt  unswept  fins,  and  has  compiled  a  fairly 
complete  bibliography  of  the  related  work. 

2.1.2  Numerical  Investigations.  The  flow  over  a  blunt  unswept  fin  appears  to 
have  been  first  solved  numerically  by  Hung,  Kordulla,  and  Bunning!46’  471  in  the  mid-1980’s 
(M=2.95,  Re/£=63x  106m-1).  They  demonstrated  good  agreement  between  their  Reynolds 
averaged  Navier-Stokes  (RANS)  solutions  and  experimental  surface  pressure  distributions 
obtained  by  Dolling  and  Bogdonoff.  t24l  Small  topology  differences  were  probably  due  to 
limited  grid  resolution. 

Soon  after,  McMaster  and  Shangl73!  numerically  simulated  this  flow  by  solving  the 
Favre  (mass)  averaged  Navier-Stokes  (FANS)  equations  with  an  algebraic  turbulence  model 
(M=2.95,  Re/^=64xl06m-1).  While  the  primary  focus  of  this  effort  was  to  investigate  the 
effect  of  fin  sweep  on  the  flowfield,  solutions  obtained  on  the  un-swept  fin  were  similar  to 
those  of  Hung  and  Bunning!46!,  with  the  separation  line  on  the  plate  agreeing  slightly  better 
with  experimental  results.  Kubendran,  et  al.!59!  obtained  similar  results,  and  designed  a 
leading-edge  fillet  that  eliminated  the  leading  edge  flow  separation. 

In  the  early  1990’s,  Knight!54'  811  reviewed  the  available  numerically-oriented  three- 
dimensional  shock/turbulent  boundary  layer  interaction  literature.  He  considered  five  basic 
interactions,  assessing  the  numerical  capability  of  the  simulation  methods  by  comparison 
with  available  experimental  data.  One  of  these  interactions  was  induced  by  a  blunt  unswept 
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fin  on  a  flat  plate.  At  the  time  of  his  article,  only  algebraic  turbulence  modeling  had  been 
employed  in  any  calculations  reported  over  this  geometry.  Also,  while  the  separation  shock 
wave  has  been  observed  to  be  an  unsteady  feature)23’  24’  251,  he  noted  that  all  computational 
simulations  of  the  supersonic  blunt  fin  problem  were  steady,  suggesting  that  this  may  have 
been  due  to  inadequate  grid  resolution  or  inaccuracies  in  the  turbulence  models.  )54J  Since 
the  experimentally  observed  unsteadiness  of  the  flowfield  involves  dominant  frequencies  on 
the  order  of  1  kHz!92!,  direct  or  large-eddy  turbulence  simulations  would  be  required  to 
resolve  the  unsteady  flowfield  computationally!46)  Only  the  time-averaged  ‘mean’  solutions 
have  been  sought  to  date,  which  are  much  more  economically  obtained. 

Since  Knight’s  review,  Chima  and  Yokota)21!,  Chen  and  Hung)29),  and  Rizzetta)88) 
among  others  have  numerically  simulated  the  flow  around  a  cylinder/flat-plate  juncture 
using  various  methods.  While  the  first  two  of  these  studies  used  algebraic  eddy  viscosity 
models  with  an  ad  hoc  modification  to  the  length  scale  to  account  for  the  multiple  inter¬ 
secting  solid  surfaces,  Rizzetta)88)  employed  the  two  equation  k-e  turbulence  model  of  Jones 
and  Launder)48) ,  including  low  Reynolds  number  terms  (M=2.5,  Re/£=19.3xl06m-1).  The 
results  agreed  very  well  with  experimental  results  -  somewhat  better  than  did  the  results 
of  Chen  and  Hung)20).  This  was  probably  due  in  part  to  higher  grid  resolution  and  strategic 
grid  clustering.  Only  one  vortex  pair  was  predicted,  while  three  pairs  were  experimentally 
observed  by  Sedney  and  Kitchens)92’ 93).  Haidinger  and  Friedrich  I37)  were  also  unable  to 
accurately  reproduce  the  upstream  effects  of  the  shock  induced  separation  using  Wilcox’s 
k-u  closure  model )118).  However,  the  downstream  predictions  of  pressure  distribution  and 
skin  friction  were  reasonable.  In  general,  recent  efforts  to  assess  the  viability  of  numerical 
methods  for  the  simulation  of  similar  complex  flows)12’  27 • 35’  49’  761  have  met  with  varying 
degrees  of  success. 

2.2  Experimental  Turbulence  Measurements 

Settles  and  Dodson)96)  recently  proposed  a  set  of  criteria  for  experimental  data  gath¬ 
ering  and  reporting.  Adherence  to  these  criteria  is  considered  necessary  in  order  for  the 
data  to  be  useful  for  CFD  research  on  supersonic  turbulent  shock/boundary  layer  interac¬ 
tions.  They  advocate  inclusion  of  each  experiment  satisfying  these  criteria  in  a  database 
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to  be  put  forth  as  a  standard  for  CFD  code  validation  and  turbulence  model  development. 
The  primary  focus  of  these  criteria  is  to  provide  to  the  CFD  community  all  of  the  necessary 
data  required  to  numerically  simulate  the  experiment.  These  criteria  have  been  adopted 
as  guidance  for  the  current  research  effort. 

As  reported  by  Settles  and  Dodson,  a  great  number  of  experiments  have  been  targeted 
at  measuring  turbulent  quantities  in  interacting  flows.  However,  most  of  these  experiments 
do  not  meet  the  high  standards  required  for  modern  code  validation.  In  fact,  at  the  date 
of  publication,  they  found  only  19  experiments  which  met  their  criteria,  none  involving 
blunt  fin-type  geometries.  Most  experiments  were  rejected  because  they  were  not  able  to 
provide  data  useful  for  testing  turbulence  models.  Others  failed  to  report  error  bounds, 
the  incoming  boundary  layer  profile,  or  simply  did  not  provide  the  data  in  a  useful  form. 
Their  survey  rejuvenated  efforts  to  obtain  such  critical  data. 

Seven  of  the  experiments  accepted  into  this  database  considered  the  interaction 
between  sharp  fin-like  geometries  and  flat  plates.  I43-  S1-  55- 60'  61  -  62> 63-  64’  65’  891  However, 
among  these  experiments,  flowfield  data  are  scarce,  and  field  turbulence  data  are  even 
more  rare.  The  current  study  represents  a  continuance  of  the  experiments  accepted  to 
date  by  introducing  surface  and  shock-generator  curvature  effects,  and  thus  provides  ex¬ 
perimental  data  that  do  not  presently  exist. 
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3.  Experimental  Methodology 


This  chapter  outlines  the  tools  and  methods  used  to  extract  meaningful  flowfield 
information.  Instruments  include  pitot  pressure  probes,  cone-static  pressure  probes,  and 
hot-film  cross-wire  probes.  The  pressure  probes  determine  mean  flow  information  such 
as  density,  pressure,  and  velocity  magnitude.  The  hot-film  cross-wire  probes  determine 
velocity  components  and  turbulence  values  in  a  form  which,  when  manipulated,  are  useful 
for  comparative  analysis  to  results  obtained  numerically. 

3.1  WAF  Model 

The  WAF  model  (Figure  3.1)  is  comprised  of  a  cylinder  of  the  fin  radius,  r=1.59cm, 
blended  to  a  removable  test  section  wall,  and  has  a  maximum  height  of  0.5r.  It  is  designed 
to  represent  a  single  fin  of  a  typical  WAF  configuration,  and  sized  to  maximize  data 
resolution  while  avoiding  tunnel  blockage.  The  fin  has  the  same  proportions  as  free-flight 
models  which  have  been  tested  at  the  Wright  Laboratory  Armament  Directorate  I1-  n3> , 
with  a  thickness  of  0.2r  (3.18mm),  a  span  of  y/2 r  (22.5mm),  and  a  chord  length,  c,  of  1.28r 
(20.3mm).  The  leading  edge  and  tip  of  the  fin  are  beveled  at  45°.  The  cylinder  is  5.12  fin 
radii  in  length  with  the  single  fin  placed  at  the  downstream  base.  Upstream,  the  cylinder 
is  blended  to  the  tunnel  floor  with  a  polynomial  chosen  to  ensure  second  order  continuity 
in  the  streamwise  direction.  The  blending  region  is  5r  long  and  starts  21.48r  from  the 
throat  of  the  wind  tunnel  nozzle. 

3.2  Mach  3  Wind  Tunnel 

The  AFIT  Mach  3  wind  tunnel  (Figure  3.2)  is  an  intermittent  blow-down  tunnel 
with  downstream  evacuation.  The  maximum  run  time  for  the  tunnel  was  30  seconds  with 
an  evacuation  time  of  6  to  10  minutes.  The  nozzle  measures  27.3cm  from  throat  to  test 
section  entrance.  The  test  section  is  33.0cm  long  with  a  square  6.35cm  cross-section.  The 
settling  chamber  total  pressure  and  temperature  are  2.0-2. 14atm  and  294K  respectively, 
yielding  a  freestream  Reynolds  number  of  Re/£=17-18xl06m-1.  Figure  3.3  is  a  photo  of 
the  nozzle  and  removable  test  sections  with  the  side-walls  removed  for  visibility. 
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Figure  3.1  Single  wrap-around  fin  model  (dimensions  in  inches). 

3.2.1  Operating  Conditions.  Detailed  pitot  and  cone- static  measurements  had 
been  previously  obtained  in  an  empty  test  section  and  used  to  determine  the  mean  inflow 
conditions  for  the  test  section.^  From  this  survey,  the  freestream  Mach  number  was  de¬ 
termined  to  be  2.9  (i7e=607m/s).  Results  from  this  earlier  experiment  were  used  to  specify 
the  upstream  boundary  conditions  for  the  preliminary  inviscid  numerical  simulation. f109! 

During  the  present  research  (with  the  model  in  the  tunnel),  the  test  section  entrance 
conditions  were  again  surveyed.  It  was  determined  that  at  a  location  just  downstream  of 
the  test  section  entrance  (8.4c  ahead  of  the  fin  leading  edge)  the  freestream  Mach  number 
was  2.80  (Ue=601m/s),  with  a  measurement  uncertainty  of  2.8%.  The  boundary  layer  on 
the  tunnel  ceiling  at  a  location  1.27cm  upstream  of  the  model  has  a  measured  thickness  of 
5.3mm  at  the  centerline,  defined  by  the  distance  from  the  surface  where  M=0.95%  Me.  A 
summary  of  measured  freestream  mean-flow  conditions  are  presented  in  Table  3.1. 
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Figure  3.3  AFIT  Mach  3  wind  tunnel. 


3.2.2  Tunnel  Instrumentation.  Instrumentation  for  the  Mach  3  tunnel  included 
an  internal  pitot  tube  and  thermocouple  in  the  plenum  chamber  to  measure  Ptoo  and  Tt<x>. 
The  thermocouple  (K-type,  Omega  part  no.  SEFE-K-5)  had  an  accuracy  of  ±lKf32h  and 
the  pitot  (total)  pressure  was  fitted  with  a  6.7atm  pressure  transducer  (Endevco  model 
851QC-1O0). 

3.2.3  Traversing  Equipment.  The  probe  traverse  system  used  in  the  Mach  3 
tunnel  used  a  stepper  motor  with  400  steps  per  revolution.  The  slide  moves  on  a  shaft  of 
40  threads  per  2.54cm,  has  a  total  travel  distance  of  16.51cm,  and  moves  at  a  maximum 
speed  of  about  0.254  cm /sec.  ^  A  photograph  of  the  traversing  assembly  is  provided  in 
Figure  3.4  (courtesy  of  Miller^74!).  Traverse  control  is  accomplished  via  personal  computer. 
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Table  3.1  Mach  3  tunnel  freestream  inflow  conditions. 


Condition 

Value 

M 

2.80 

Pt 

217kPa  =  2.14atm 

P 

7.99kPa  =  0.0789atm 

Tt 

294K 

f 

114.5K 

u 

p 

0-2431^: 

Figure  3.4  Traverse  assembly  (Miller). 
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Figure  3.5  AFIT  Mark  5  wind  tunnel  with  WAF  model. 


3.3  Mach  5  Wind  Tunnel 

Experimental  results  were  also  obtained  in  the  new  AFIT  Mach  5  wind  tunnel  (Fig¬ 
ure  D.l).  This  blow-down  wind  tunnel  has  a  heated  air  supply  and  has  been  operated  over 
a  range  of  unit  Reynolds  numbers  (Re/^«32-75xl06m-1).  Although  the  tunnel  is  capable 
of  producing  much  longer  run-times,  the  current  experiments  only  required  total  run-times 
of  10-15  seconds.  The  design,  construction,  instrumentation,  and  initial  calibration  of  this 
facility  comprised  a  significant  portion  of  this  research  effort.  More  detailed  information 
on  the  AFIT  Mach  5  wind  tunnel  is  available  in  Appendix  D. 

S.Jf  Pressure  Probes 

A  pitot  probe  and  a  10°  (±.03°)  cone-static  probe  were  used  to  measure  pressures 
throughout  the  flowfield  in  the  Mach  3  wind  tunnel  (Figure  3.6).  Total  and  static  pres¬ 
sure  data  are  obtained  directly  from  the  pitot  and  cone-static  probes  through  normal  and 
axisymmetric  oblique  shock  relations.  These  data  can  be  used  to  calculate  mean  flow 

information  such  as  Mach  number,  density,  and  the  magnitude  of  the  velocity.  Through 

- 

. 
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Figure  3.6  Pressure  probes. 


the  manipulation  of  results  from  conical f9’  56i  and  normal  shock  relations!10!,  the  following 
curve  fit  was  proposed  by  Bowersox.E1^ 

&  =  -0.052976  +  4.6840fio°  -  18.6786^20o  +  50.7006&C  -  54.1577^  (3-1) 

where  fio<>  is  defined  as  the  ratio  of  the  pressure  measured  with  the  10°  cone-static  probe 
to  the  pitot  pressure  (6o°=Fc/Pt2=/cn[M]).  It  has  been  found  that  for  flowr  angles  less 
than  6.0°  the  errors  in  Mach  number  axe  less  than  0.03J113!  Also  note  Equation  3.1  is  valid 
for  Mach  numbers  in  the  range  from  1.5  to  4.4  and  has  a  standard  deviation  of  0.06%. I13! 

Stronger  pitot  and  cone-static  probes  than  those  used  in  the  Mach  3  tunnel  were 
fabricated  from  heavy  gauge  stainless  steel  to  withstand  the  higher  loads  associated  with 
the  Mach  5  tunnel.  This  cone-static  probe  had  a  semi-vertex  angle  of  20°,  which  produces 
higher  measured  pressure  levels  than  the  10°  probe,  enabling  more  accurate  readings.  The 
higher  cone  angle  also  minimized  the  error  due  to  flow’  angularity.  For  this  new’  cone  angle, 
and  to  envelope  the  higher  Mach  number  range,  a  new  relationship  between  Mach  number 

and  the  measured  pressures  was  needed.  Again,  conical  and  normal  shock  solutions  were 

• 

used  to  develop  a  direct  relationship  between  Mach  number  and  the  ratio  of  measured 
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pressures. 

fa  =  0.3469176  +  0.908347^20°  -  0.0518062/^0°  (3-2) 

This  relationship  is  valid  for  Mach  numbers  between  1.5  and  6.5  with  a  correlation  co¬ 
efficient  of  99.990%.  Note  that  this  form  provides  monotonic  behavior  at  high  pressure 
ratios  (low  Mach  numbers),  and  is  nearly  as  accurate  within  the  range  of  validity  as  a 
fourth-order  polynomial. 

With  the  exception  of  the  Mach  5  pitot  pressure,  all  probes  were  connected  with 
tubing  to  a  latm  transducer  (Endevco  model  8510C-15),  which  sent  its  signal  to  a  signal 
processor  (Endevco  model  4428A).  This  unit  uses  the  maximum  range  of  the  pressure  trans¬ 
ducer,  is  self-zeroing,  and  connects  directly  to  the  data  acquisition  system  (Section  3.6). 
The  Mach  5  pitot  pressure  was  sensed  with  a  6.7atm  transducer  (Endevco  model  8510C- 
100).  Atmospheric  pressure  was  recorded  using  a  Druck  resonant  sensor  barometer. 

3.5  Hot- Film  Cross- Wire  Probes 

Presented  here  is  a  brief  overview  of  the  hot-film  cross- wire  probes  and  the  associated 
data  reduction  techniques.  More  details  on  the  hot-film  anemometry  methods  employed 
to  determine  both  mean  and  turbulent  flow  information  are  provided  in  Appendix  B. 

Two  hot-film  probes  (Figure  3.7)  were  used  to  determine  the  velocity  and  fluctuations 
of  the  flow  in  the  Mach  3  tunnel.  Both  probes  were  two-component  cross-wire  hot-film 
probes,  each  with  two  thin  films  of  platinum  1mm  long  and  51^m  in  diameter.  One  probe 
had  the  films  oriented  in  the  vertical  plane,  angled  at  ±45°  to  the  horizontal  (u-v)  plane. 
The  other  probe’s  films  were  oriented  in  the  horizontal  plane  and  similarly  angled  (u-w 
plane).  The  transverse  separation  between  the  two  films  was  1mm.  During  most  runs, 
these  probes  were  traversed  about  5cm.  Since  the  traverse  was  moving  slowly  relative 
to  the  flow,  pressure  and  hot-film  measurements  were  taken  while  the  traverse  was  in 
motion.  This  technique  has  been  validated  previously  in  several  experiments  in  the  Mach  3 
tunnel.  t26> 69'  74  ■  75 ■  941  The  hot-film  probes  were  connected  to  a  TSI  brand  Intelligent  Flow 
Analyzer!111].  A  Tektronix  2454B  oscilloscope  was  used  to  view  and  tune  the  frequency 
response  shape  and  response  time  for  each  film  on  the  probe. 
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[a]  u-v  TSI  model  1243-20. 


[b]  u-w  TSI  model  1243AN-20. 
Figure  3.7  Hot  film  probes. 


Cross-film  measurements  provide  a  means  to  resolve  the  mean  mass  flux  vector  into 
its  Cartesian  components.  When  the  Mach  number  normal  to  the  film  is  greater  than  1.2, 
the  Prandtl  number  may  be  assumed  constant!57'  981.  When  the  aspect  ratio  of  the  hot-film 
is  also  much  larger  than  unity,  then  the  functional  form  of  the  Nusselt  number  may  be 
further  simplified.  It  has  been  determined  experimentally  that  King’s  Law ,  the  functional 
relationship  between  Nusselt  number  and  effective  Reynolds  number  for  incompressible 
flow,  is  also  an  acceptable  relationship  for  compressible  flow.!14'  13>  161 

Nu  =  <z*\/Re7+  bk  (3.3) 

The  unknowns  in  Equation  3.3  are  the  calibration  coefficients,  a *  and  and  Ree  which 
is  the  effective  cooling  Reynolds  number  normal  to  the  film  based  on  film  diameter  and 
reference  viscosity.  The  hot-film  calibration  is  performed  by  placing  the  probe  in  the  tunnel 
free  stream,  varying  the  plenum  pressure  (thus  changing  Ree),  and  measuring  the  hot-film 
voltage,  providing  Nu.  The  calibration  constants  a*,  and  bk  are  then  determined  using  a 
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least  squares  linear  regression  on  the  data.  Hence,  the  calibrated  hot-film  provides  values 
of  local  Reynolds  number,  essentially  allowing  direct  measurement  of  the  mass-flux,  pU . 

A  complication  faced  by  modelers  is  that  the  apparent  mass  terms  appearing  in  the 
compressible  turbulence  terms  of  the  RANS  equations  (Appendix  A)  cannot  be  directly 
measured  with  conventional  methods.  One  of  the  key  features  of  cross-wire  anemometry 
is  how  the  Reynolds  turbulent  shear  stress  can  be  directly  measured  as  the  negative  of  the 
mass-flux  correlation  term  combined  with  the  density  fluctuation.  Bowersox  and  Schetzl16! 
have  shown  that  these  apparent  mass  terms  may  be  approximated  by  first  noting  that 

p<t>  =  p<p  +  ( p4>)'  (3.4) 


and  also 


p<t>  =  (p  +  p')(4>  +  4>)  =  p<t>  +  p4>  +  p'<t>  +  p'4>  (3.5) 


Then,  using  (3.4)  and  (3.5)  in  the  Navier-Stokes  equations  (A.8)  and  neglecting  third  order 
terms  allows  us  to  rewrite  the  Reynolds  shear  stress  (Equation  A.ll)  as 


(gj 


(3.6) 


where  the  second  term  has  been  shown  to  be  much  less  than  the  first  term  for  thin  layer 
flows.  In  the  present  experiments,  the  second  term  was  always  at  least  an  order  of  magni¬ 
tude  smaller  than  the  first  term.  The  first  term  may  be  directly  measured  with  cross-film 
probes. 

Multiple  overheat  ratios  are  generally  required  to  determine  all  of  the  flow  informa¬ 
tion.  However,  sensitivity  to  total  temperature  fluctuation  is  often  negligible  for  hot-films 
operating  at  high  overheat  ratios  (j^  >  2.0).  In  cases  where  the  fluctuation  in  total 
temperature  is  negligible,  only  a  single  overheat  ratio  is  necessary  to  determine  flow  field 
information.  Since  the  AFIT  Mach  3  wind  tunnel  has  been  found  to  maintain  total  temper¬ 
ature  fluctuation  below  about  2.0%.  I44- 45-  70 ■  751,  a  single-overheat  data  reduction  method 
was  used.  The  single-overheat  hot-film  analysis  and  its  use  in  the  present  research  is 
detailed  in  Appendix  B. 
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3.6  Data  Acquisition 

The  AFIT  Mach  3  tunnel  is  instrumented  with  a  Multiprot79!  data  acquisition  system. 
The  Multipro  has  four  Model  120  data  acquisition  boards  installed  with  one  megabyte  of 
memory  per  board.  The  selection  of  the  sampling  frequency  was  based  in  part  on  the 
criteria  of  Smits  and  Mucld97! .  This  sampling  frequency  was  selected  such  that  the  character 
of  the  flow  was  captured  (as  indicated  by  the  sampling  energy  spectra).  Experiments 
were  conducted  to  insure  that  data  were  recorded  at  a  sufficiently  high  frequency.  Data 
acquisition  parameters  are  summarized  in  Table  3.2. 


Table  3.2  Data  acquisition  system  parameters. 


Acquisition  Parameter 

M=2.8  Experiment 

M=4.9  Experiment 

Pressure  Probe  Sampling  Frequency 

200Hz 

500Hz 

Traverse  Speed 

.077in/sec 

.20in/sec 

Hot-film  Response  Time 

8-12  fis 

— 

Hot-film  Sampling  Frequency 

16.6kHz 

— 

3. 7  Shadowgraph  and  Schlieren  Optics 

The  shadowgraph  and  schlieren  images  were  obtained  on  Polaroid  film  using  a  Xenont121) 
arc  light  source  with  a  spark  duration  of  10  nanoseconds.  For  the  shadowgraphs,  the  light 
is  reflected  from  a  101.6  cm  focal  length  mirror  and  photographed  with  series  52  Polaroid 
film.  The  schlieren  imaging  requires  an  additional  mirror  and  knife  edge  to  polarize  the 
incoming  light  to  the  film.  The  diagrams  for  the  shadowgraph  and  schlieren  layout  are 
provided  in  Figures  3.8  and  3.9. 
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Figure  3.9  Schlieren  camera  and  mirror  setup. 
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4-  Numerical  Methodology 


In  this  chapter,  the  numerical  methods  which  have  been  applied  in  the  present  re¬ 
search  are  presented,  along  with  the  associated  boundary  conditions.  Solutions  to  the 
governing  equations  (Appendix  A)  were  obtained  using  the  General  Aerodynamic  Sim¬ 
ulation  Program  (GASP).  !7’  72>  GASP,  a  fully  conservative  shock  capturing  code,  has 
been  widely  used  by  the  CFD  community  for  the  analysis  of  supersonic  and  hypersonic 
flows.!36'  71  - 90’  1141  GASP  Version  2.0  was  used  for  the  inviscid  calculations,  and  GASP 
Version  3.06  was  used  for  the  viscous  calculations.  A  sample  GASP  3.0  input  deck  for  the 
viscous  simulation  at  Mach  2.8  in  included  in  Appendix  E. 

4.1  Inviscid  Calculations  (M=2.9,5.0) 

As  a  means  of  establishing  the  most  suitable  locations  for  experimental  measure¬ 
ments,  pilot  numerical  simulations  were  conducted  by  solving  the  Euler  equations  upstream 
and  in  the  vicinity  of  the  WAF.!109)  These  simulations  were  conducted  on  grids  representing 
the  WAF  model  as  installed  in  the  Mach  3  and  Mach  5  tunnel  test  sections,  and  will  be 
discussed  further  here. 

4.1.1  Grid  Definition.  To  resolve  the  features  of  the  shock  structure  in  the  vicin¬ 
ity  of  the  fin,  the  grid  was  clustered  near  the  fin  and  cylinder  surfaces.  For  both  Mach 
numbers,  the  entire  computational  mesh  consists  of  11  computational  zones  connected  via 
19  zonal  boundaries.  Flow  variable  values  are  passed  through  the  zonal  boundaries  via 
a  simple  five-point  overlap.  The  entire  test  section  (all  four  walls)  was  included  in  the 
Mach  2.9  simulation.  Since  at  Mach  5.0  the  shock  emanating  from  the  body  was  expected 
to  remain  very  close  to  the  fin,  only  two-thirds  of  the  test  section  was  modeled.  This  sim¬ 
plification  provided  substantial  savings  in  computational  time  and  memory  requirements. 

The  faces  of  each  zone  in  the  region  near  the  fin  were  defined  with  an  elliptic  grid 
generator  (GRIDGENI85))  employing  Thomas  and  Middlecoff!106)  control  functions.  Grid 
orthogonality  was  enforced  along  the  edges  of  each  zone,  while  constraining  the  grid  points 
to  lie  on  the  geometry  surface.!101- 100]  The  interior  points  of  each  zone  were  then  positioned 
using  a  standard  trans-finite  interpolation  (TFI)  scheme.!85)  Figure  4.1  shows  the  numerical 
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Figure  4.1  Grid  boundaries  and  zonal  structure  for  inviscid  Mach  2.9  simulation. 


representation  of  the  Mach  3  tunnel  test  section  near  the  fin.  The  Mach  5.0  simulation 
was  performed  on  a  grid  created  by  the  same  methods  and  defined  by  the  same  model 
geometry.  However,  the  tunnel  side  walls  were  moved  to  represent  the  wider  test  section. 

4*1.2  Solution  Strategy.  Since  the  flow  disturbances  do  not  propagate  upstream 
in  a  fully  supersonic  flow,  the  solution  was  space-marched  from  an  upstream  starting  plane 
0.8r  ahead  of  the  blended  body  to  a  location  approximately  0.4c  ahead  of  the  fin  (see 
Figure  5.1).  The  location  at  which  the  space  marching  was  concluded  was  chosen  to  be 
well  ahead  of  the  expected  bow  shock  location.  Pilot  2-D  simulations  over  a  beveled  fin 
provided  early  estimates  of  the  position  of  the  shock  induced  by  the  fin.  These  simulations 
also  provided  an  expedient  means  to  ensure  that  the  grid  density  and  clustering  near  the 
fin  to  be  used  in  the  3-D  simulation  would  be  sufficient  to  capture  the  key  feature  of  the 
detached  fin  shock. 


Table  4.1  Freestream  conditions  used  for  Euler  calculations. 


Condition 

Value 

Value 

M 

2.9 

5.0 

Pt 

20.4atm 

P 

0.0633atm 

0.03858atm 

Tt 

294K 

375K 

T 

109.6K 

62.5K 

U 

■si 

P 

■mu 

■mon 

The  upstream  boundary  condition  (initial  marching  condition)  was  simply  prescribed 
to  be  the  nominal  Mach  3  wind  tunnel  freestream  condition  (751  summarized  in  Table  4.1. 
The  upstream  boundary  condition  for  the  Mach  5.0  simulation  was  set  at  the  tunnel 
freestream  design  conditions  (see  Appendix  D).  Flow  tangency  was  prescribed  on  all  solid 
surfaces,  and  was  enforced  explicitly  to  first-order  accuracy  with  a  full  flux  method.  t72) 
For  the  Mach  5.0  case,  flow  conditions  were  extrapolated  (first-order)  from  the  interior  at 
the  upper  boundary  of  the  computational  domain  (opposite  the  model).  The  two-factor 
approximately-factorized  equations  were  solved  at  each  streamwise  marching  plane,  with 
relaxation  applied  in  the  marching  direction.  A  third  order  upwind  biasing  (full  flux) 
scheme  was  employed  in  the  marching  direction,  while  a  second  order  up-winding  scheme 
with  Roe’s  flux  difference  splitting  was  used  in  the  cross-flow  plane. 

Beginning  at  a  location  0.4c  ahead  of  the  fin,  the  three-factor  approximately-factorized 
equations  were  solved  globally  to  first-order  spatial  accuracy  by  employing  van  Leer  split¬ 
ting  of  the  inviscid  fluxes!112! .  The  upstream  boundary  condition  for  this  region  is  specified 
to  be  the  solution  given  at  that  plane  by  space-marching.  Flow  conditions  were  extrapo¬ 
lated  to  first-order  from  the  interior  at  the  exit  boundary  of  the  computational  domain, 
which  is  at  the  trailing  edge  of  the  WAF  model.  Again,  flow  tangency  is  enforced  on  all 
solid  surfaces  (including  the  tunnel  walls)  explicitly  to  first-order  accuracy  with  a  full  flux 
method.  For  the  Mach  2.9  case,  this  region  is  comprised  of  8  computational  zones,  con¬ 
taining  a  total  of  825xl03  cells.  For  the  Mach  5.0  case,  this  region  is  also  comprised  of  8 
computational  zones,  containing  a  total  of  655 xlO3  cells. 
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4-1-3  Convergence  Issues.  The  inviscid  calculations  were  not  intended  to  be  a 
rigorous  computational  study.  The  primary  goal  of  the  inviscid  study  was  to  obtain  a 
preliminary  estimate  of  acceptable  instrument  locations  for  the  experimental  study.  The 
inviscid  computations  were  assumed  to  be  converged  when  the  residual  of  the  state  vector 
was  reduced  by  three  orders  of  magnitude.  At  a  CFL  number  of  1.0®,  this  required  838 
iterations  for  the  Mach  2.9  calculation  and  643  iterations  for  the  Mach  5.0  calculation. 

4-1-4  Computational  Requirements.  The  inviscid  computations  were  performed 
on  a  Cray  C916/161024  supercomputer.  The  global  calculations  near  the  fin  required  less 
that  16mWords  of  memory  and  23.5  seconds  per  iteration  to  complete.  The  calculations 
on  the  slightly  smaller  Mach  5.0  grid  required  17.5  seconds  per  iteration. 

4-1-5  Inviscid  Fin  Simulations  at  Other  Mach  Numbers.  The  flow  near  the  fin  has 
also  been  numerically  investigated  at  several  Mach  numbers  other  than  2.9  and  5.0  using 
the  same  methods  discussed  above.  For  these  simulations,  the  effect  of  the  blending  region 
and  tunnel  were  neglected  by  imposing  the  freestream  boundary  condition  shortly  upstream 
of  the  fin.  This  assumption  minimized  computational  and  solution  storage  requirements. 
For  this  study,  the  computational  grid  was  identical  to  the  8-zones  near  the  fin  from  the 
inviscid  Mach  2.9  study  (Figure  4.1). 

4-2  Viscous  Calculations 

The  Reynolds-averaged  Navier-Stokes  (RANS)  equations  were  solved  over  the  model 
geometry  using  GASP  version  S.OSt7!,  with  the  algebraic  turbulence  model  of  Baldwin  and 
Lomax  I11). 

4-2.1  Grid  Definition.  The  numerical  representation  of  the  test  article  used  for 
the  viscous  calculations  is  provided  in  Figure  4.2.  A  multi-zone  approach  was  again  taken. 
The  entire  computational  mesh  consisted  of  12  computational  zones  (the  edges  of  which 
are  shown  in  Figure  4.2)  connected  by  21  zonal  boundaries  and  was  comprised  of  8.2 xlO6 
cells.  The  flow  variable  values  were  passed  through  the  zonal  boundaries  via  a  simple 
five-point  overlap. 
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extrapolation 


Figure  4.2  Grid  boundaries  and  zonal  structure  for  viscous  simulations. 

To  resolve  the  features  of  the  flowfield  in  the  vicinity  of  the  fin  and  to  provide  the 
resolution  required  by  turbulence  models,  the  grid  was  clustered  near  the  fin  and  cylinder 
surfaces.  At  a  location  0.4c  upstream  of  the  fin  leading  edge  the  wall  grid  spacing  corre¬ 
sponded  to  a  y+  value  of  roughly  0.15  for  both  Mach  numbers  investigated.  At  this  location, 
approximately  half  of  the  points  in  the  normal  direction  were  contained  in  the  boundary 
layer.  Zone  faces  near  the  fin  were  defined  using  an  elliptic  grid  generator  (GRIDGENl85!) 
employing  Thomas  and  MiddlecofP106]  control  functions.  Grid  orthogonality  was  enforced 
along  the  edges  of  each  zone,  while  constraining  the  grid  points  to  lie  on  the  geometry 
surface.  I101-  10°l  The  interior  points  of  each  zone  were  then  positioned  using  a  standard 
trans-finite  interpolation  (TFI)  scheme.  I85) 
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4-2.2  Solution  Strategy.  For  all  calculations  the  thin  layer  viscous  terms  (with 
Sutherland’s  viscosity  model)  were  included  in  the  governing  equations.  The  algebraic  eddy 
viscosity  model  of  Baldwin  and  Lomax!11!  was  also  used  to  simulate  turbulence.  Details  of 
this  turbulence  model  are  available  in  A. 5.  Since  flow  disturbances  do  not  propagate  up¬ 
stream  in  a  fully  supersonic  flow,  and  since  experimental  results  indicated  that  the  blended 
body  produced  no  separated  flow  regions,  the  parabolized  Navier- Stokes  (PNS)  equations 
were  solved  to  a  location  0.5c  ahead  of  the  fin.  This  location  was  deemed  sufficiently  far 
upstream  of  the  fin  interaction  region  based  on  previously  conducted  visualization  experi¬ 
ments.!45-  109] 

To  allow  for  the  specification  of  a  two-dimensional  upstream  boundary  condition,  the 
PNS  equations  were  solved  on  a  two-dimensional  grid  for  a  short  distance  (0.8r)  ahead  of 
the  model  body.  At  the  inflow  boundary,  an  experimental  profile  determined  using  pitot, 
cone-static,  and  hot-film  probes  was  prescribed  (see  Appendix  E).  Owing  to  probe  volume 
effects,  experimental  data  could  not  be  obtained  sufficiently  close  to  the  wall  to  include 
the  laminar  sublayer.  However,  the  two-dimensional  PNS  region  allowed  the  boundary 
layer  to  develop  into  a  fully  turbulent  profile  upstream  of  the  blended  body  region.  The 
two-dimensionality  of  the  flow  in  the  AFIT  Mach  3  wind  tunnel  has  been  well  docu¬ 
mented.  f75-  109l  For  the  simulation  of  the  fin  in  the  Mach  5  wind  tunnel,  inflow  conditions 
were  derived  from  pitot  pressure  data.  The  two- dimensionality  of  the  flow  in  the  AFIT 
Mach  5  tunnel  is  documented  in  Appendix  D. 

The  two-dimensional  solution  was  mapped  to  the  three-dimensional  grid  at  the  lead¬ 
ing  edge  of  the  blended  body.  As  in  the  inviscid  simulations,  the  two- factor  approximately- 
factorized  equations  were  solved  at  each  marching  plane  by  employing  a  third  order  upwind 
biasing  scheme  with  relaxation  in  the  marching  direction  and  a  second  order  up-winding 
scheme  with  Roe’s  flux  difference  splitting  in  the  cross-flow  plane.  The  symmetry  of  the 
model  body  was  exploited  by  solving  the  PNS  equations  over  only  half  of  the  blended 
body  region  and  by  employing  an  x-y  symmetry  condition  at  the  z=0  plane.  GASP  was 
modified  to  allow  the  solution  at  the  final  symmetric  marching  plane  to  be  reflected  across 
x-y  plane.  The  solution  at  this  plane  was  then  used  as  the  upstream  condition  for  the 
downstream  asymmetric  region.  In  the  vicinity  of  the  fin  (behind  the  plane  at  0.5c  ahead  of 
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the  fin)  the  RANS  equations  were  solved  to  third-order  spatial  accuracy  using  Jacobi  inner 
iterations!7!.  The  inviscid  fluxes  were  split  by  the  method  of  van  Leer!112!,  and  the  min-mod 
limiter!40!  was  used.  This  region  is  comprised  of  8  computational  zones,  containing  a  total 
of  4.2 xlO6  cells. 

On  the  model  surfaces,  shown  as  a  mesh  in  Figure  4.2,  a  no-slip  condition  on  the  ve¬ 
locity,  an  isothermal  wall  temperature,  and  vanishing  normal  pressure  gradient  (boundary 
layer  approximation)  were  enforced. 

dP 

iij  =  0  T  —  T-wall  ~dn  —  ^ 

The  wall  temperature  was  294K  for  the  Mach  2.8  simulation.  For  the  Mach  4.9  simulation, 
the  effect  of  wall  temperature  was  investigated  by  enforcing  relatively  cold  (294K)  and 
hot  (340K)  isothermal  boundary  conditions  at  the  wall.  An  adiabatic  wall  condition  was 
also  examined.  From  these  simulations,  it  was  determined  that  the  thermal  boundary 
condition  at  the  wall  has  a  negligible  influence  on  the  solution.  The  results  presented  here 
were  obtained  with  a  294K  isothermal  boundary  condition,  and  are  indistinguishable  from 
those  obtained  using  the  other  thermal  boundary  conditions. 

The  solid  surface  boundary  conditions  were  enforced  explicitly  to  second-order  ac¬ 
curacy  with  a  full  flux  method!7!.  Based  on  the  results  of  the  previously  conducted  ex¬ 
perimental  and  inviscid  numerical  investigations!45' 109! ,  it  was  presumed  that  the  side  and 
opposing  tunnel  walls  had  a  minimal  influence  on  the  flowfield  near  the  fin.  Thus,  al¬ 
though  the  size  of  the  computational  domain  represents  the  test  section,  flow  conditions 
were  extrapolated  (first-order)  from  the  interior  at  these  boundaries.  This  afforded  great 
computational  savings  by  avoiding  the  need  to  enforce  a  no-slip  condition  on  the  opposite 
and  side  walls.  Also,  since  the  side  and  opposing  tunnel  walls  were  not  modeled,  the  same 
grid  was  used  for  both  the  Mach  2.8  and  Mach  4.9  calculations.  The  flow  conditions  were 
also  extrapolated  at  the  downstream  plane. 

A  substantial  savings  of  computational  resources  was  realized  through  the  use  of  a 
grid  sequencing  method.  This  method  reduces  the  mesh  resolution  by  combining  every  N 
computational  cell  edges  into  one  cell  edge  in  a  given  direction,  where  N  is  the  sequencing 
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level.  Therefore  the  cellular  dimension  of  the  finer  grid  in  the  sequenced  direction  must 
be  divisible  by  N.  This  coarser  grid  may  again  be  sequenced  if  it  has  suitable  dimensions. 
In  the  present  research,  the  finest  grid  was  dimensioned  such  that  it  could  be  sequenced 
twice  with  N=2  in  all  three  grid  directions.  This  simply  means  that  the  dimension  of 
every  computational  block  and  of  every  face  which  defined  a  geometry  had  to  be  evenly 
divisible  by  both  2  and  4.  Initial  iterations  were  performed  on  the  ‘coarse’  grid,  then  that 
solution  interpolated  to  the  ‘medium’  and  the  ‘fine’  grids.  This  resulted  in  a  huge  savings 
of  computational  effort,  as  each  iteration  on  a  N= 2  sequenced  grid  requires  only  about 
one-eighth  the  CPU  time  and  memory  as  an  iteration  on  its  parent  grid.  Comparisons 
between  the  solution  on  the  sequenced  grids  were  made  to  investigate  grid  consistency  and 
convergence  issues. 

4-2.3  Convergence  Issues. 

4-2.3. 1  Grid  Consistency.  Grid  sequencing  not  only  accelerates  solution 
convergence  in  almost  all  temporally  integrated  problems,  but  it  also  affords  an  expedient 
means  to  evaluate  grid  consistency.  Grid  sequencing  was  also  employed  in  the  upstream 
region  which  was  evaluated  by  space-marching  the  PNS  equations  (x<— 0.5c)  as  a  means 
to  determine  grid  consistency.  Solutions  indicate  (Figure  4.3)  that  the  boundary  layer 
predicted  on  the  ‘medium’  grid  on  centerline  at  the  exit  plane  of  the  space  marched  region 
is  unchanged  by  further  grid  refinement.  Data  from  the  fine-grid  solution  was  used  to 
establish  the  upstream  boundary  condition  for  the  adjoining  region  which  surrounded  the 
fin. 

Although  the  grid  was  defined  with  sufficient  resolution  at  the  solid  surfaces  to  em¬ 
ploy  the  algebraic  turbulence  model  on  the  finest  grid,  the  wall  spacing  increases  with 
sequencing,  possibly  leading  to  erroneous  eddy  viscosity  values.  Since  the  implementation 
of  algebraic  closure  models  in  regions  having  more  than  one  physical  length  scale  (i.e.  near 
two  or  more  solid  surfaces)  may  cause  convergence  difficulties,  convergence  studies  were 
conducted  by  solving  the  laminar  form  of  the  RANS  equations  in  the  fin  region. 
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Figure  4.3  Effect  of  grid  refinement  on  predicted  boundary  layer  profile  at  the  exit  of  the 
space-marched  region  (z=0,®=— .5c). 

Pitot  pressures  from  converged  solutions  on  each  sequence  are  presented  in  Figure  4.4 
with  experimental  data  obtained  at  locations  aft  of  the  bow  shock.  The  solutions  provided 
by  the  two  finest  grids  both  exhibit  the  same  physical  characteristics  such  as  the  pitot 
pressure  increase  in  the  boundary  layer  on  the  convex  side  of  the  fin  (Y /6o<0.5)  and  a 
small  separation  near  the  body  on  the  concave  side  of  the  fin.  On  the  convex  side  of  the 
fin,  the  predicted  pitot  pressures  become  more  similar  to  the  experimental  data  with  grid 
refinement.  On  the  concave  side  of  the  fin,  the  pressure  profiles  computed  on  the  two  finest 
grids  are  very  similar,  but  the  extent  of  the  flattening  of  the  pitot  pressures  near  the  wall  is 
underpredicted,  and  the  pitot  pressure  levels  in  the  inner  boundary  layer  are  overpredicted. 
As  will  be  shown  in  Chapter  5,  the  extent  of  the  separation  region  on  the  concave  side  of 
the  fin  (Figure  4.4[b])  was  predicted  more  accurately  with  the  subsequent  addition  of  an 
algebraic  turbulence  model.  It  will  also  be  shown  that  there  are  very  large  turbulent  shear 
stresses  on  the  concave  side  of  the  fin  are  probably  responsible  for  the  overprediction  of 
the  momentum  levels  in  the  boundary  layer. 

Examination  of  the  vortical  structures  ahead  of  the  fin  leading  edge  (Figure  4.5) 
indicates  that  the  structure  of  the  flowfield  obtained  on  the  coarsest  grid  is  different  than 
that  predicted  on  the  two  finer  grids.  The  first  grid  refinement  provides  enough  resolution 
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Figure  4.4  Effect  of  grid  refinement  on  pitot  pressure  at  downstream  location. 

to  predict  a  4-vortex  structure.  Further  grid  refinement  moves  the  predicted  locations  of 
these  vortices  slightly  toward  the  fin.  It  is  important  to  note  that  this  region  is  the  most 
difficult  area  of  the  flowfield  to  resolve,  since  the  flow  character  is  changing  drastically 
over  short  distances.  However,  previous  numerical  studies  on  blunt  fins  and  cylinders  have 
indicated  that  it  may  not  be  necessary  to  resolve  this  region  of  the  flow  precisely.  This 
is  primarily  due  to  the  fact  that  these  particular  vortical  structures  are  swept  away  from 
the  fin,  minimizing  the  effects  of  their  inner  composition  on  the  fin  loading.  Indeed,  the 
agreement  can  be  very  good  between  calculated  and  experimental  quantities  away  from 
this  region  such  as  surface  pressure  on  the  fin  and  on  the  body  near  the  fin,  even  when  the 
number  of  vortices  within  the  structure  is  incorrectly  predicted. t37’  881  In  the  present  case, 
while  the  pressures  gradients  ahead  of  the  fin  on  the  model  surface  are  slightly  elevated  by 
the  last  grid  refinement,  the  pressures  on  the  fin  surface  were  virtually  identical.  Predicted 
rolling  moment  coefficients  calculated  by  integrating  the  surface  pressures  on  the  fin  sides 
for  the  two  finest  meshes  are  0.002077  (medium)  and  0.002074  (fine),  indicating  that  grid 
convergence  was  achieved  in  the  region  of  interest.  This  moment  coefficient  is  about  the 
virtual  missile  x-axis  and  based  on  body  diameter  and  cross-sectional  area. 


4-10 


xfc  xfc  xfc 

[a]  Coarse  grid.  [b]  Medium  grid.  fc]  Fine  grid. 

Figure  4.5  Effect  of  grid  resolution  on  the  computed  vortical  structure  ahead  of  the  fin 
leading  edge:  M=2.8  laminar  RANS  simulation. 

4. 2.3.2  Temporal  Convergence.  For  the  laminar  calculations,  the  two  coars¬ 
est  sequences  were  time-integrated  until  R^clO-4,  where  Rj^  is  defined  as  the  Euclidean 
(£2)  norm  of  the  residual  vector,  and  is  normalized  its  value  on  the  first  iteration  of  the 
current  sequence  (Figure  4.6).  When  approaching  convergence  on  the  finest  grid,  the  resid¬ 
ual  vector  was  dominated  by  fluctuations  in  the  vicinity  of  shock  waves.  This  behavior  is 
exhibited  in  many  problems  where  the  fiowfield  is  comprised  of  regions  in  which  the  fluid 
motion  occurs  on  greatly  disparate  time  scales.  In  such  simulations,  convergence  must,  be 
based  on  some  physical  properties  of  the  flow.  For  the  current  problem,  the  pitot  pressure 
on  lines  corresponding  to  the  experimental  data  was  chosen  to  determine  temporal  conver¬ 
gence.  The  turbulent  calculations  were  performed  on  the  finest  grid  only  with  the  laminar 
solution  as  the  Initial  condition.  Temporal  convergence  was  established  in  the  same  fashion 
as  it  was  for  the  laminar  calculation. 

4-2-4  Mach  4-9  Calculations.  The  Mach  4.9  calculation  was  conducted  following 
the  same  strategy  used  for  the  Mach  2.8  investigation,  except  that  no  laminar  computations 
were  performed.  The  Baldwin-Lomax  turbulence  model  was  used  throughout  the  simula¬ 
tion,  including  the  solutions  on  the  sequenced  grids.  The  convergence  behavior  observed 
during  these  calculations  was  very  similar  to  that  observed  in  the  laminar  computations 
at  Mach  2.8  (Figure  4.7).  Changes  in  the  vortical  flow  structure  ahead  of  the  fin  due 
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Figure  4.6  Residual  history  for  M=2.8  laminar  calculations  in  global  region. 


Iteration 

Figure  4.7  Residual  history  for  turbulent  M=4.9  calculations  in  global  region. 

to  grid  refinement  (Figure  4.8)  were  also  similar  to  those  seen  in  the  laminar  Mach  2.8 
computations  (Figure  4.5). 

4-2.5  Computational  Requirements.  The  computational  requirements  for  all  of 
the  RANS  calculations  are  provided  in  Tables  4.2  and  4.3.  As  expected,  the  computational 
requirements,  both  in  terms  of  memory  and  CPU,  increase  by  approximately  a  factor  of 
eight  with  each  grid  refinement.  The  Baldwin-Lomax  turbulence  model  required  about 
5.5%  extra  memory,  and  34%  more  time  per  iteration.  All  viscous  computations  were 
performed  on  a  Cray  C916  supercomputer. 
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Figure  4.8  Effect  of  grid  resolution  on  the  computed  vortical  structure  ahead  of  the  fin 
leading  edge:  M=4.9  RANS  simulation  with  Baldwin-Lomax. 


Table  4.2  Computational  requirements  for  Mach  2.8  RANS  calculations. 


Solution 

Method 


Laminar 
Laminar 
Laminar 
B-L  turb 


Grid 

Sequence 


3 

2 

1 

1 


Number  Required 
of  Cells  Memory 


65168  3.9mW 

521344  15.3mW 

4170752  :  107.mW 
4170752  113.mW 


Table  4.3  Computational  requirements  for  Mach  4.9  RANS  calculations. 


Solution 

Method 

Grid 

Sequence 

Number 
of  Cells 

Required  j 
Memory  j 

Time 

(sec/itn) 

B-L  turb 
B-L  turb 
B-L  turb 


65168 

521344 

4170752 


4.1mW 

16.1mW 

113-mW 


5.  WAF  Flowfield  at  Mach  2.8 


The  mean  and  turbulent  structure  of  the  flowfield  around  a  single  wrap-around  fin 
mounted  on  a  semi-cylindrical  body  has  been  investigated  using  the  experimental  and 
numerical  methods  outlined  in  Chapters  3  and  4.  In  this  chapter,  the  results  of  the  inves¬ 
tigations  of  the  WAF  at  Mach  2.8,  including  an  analysis  of  the  flowfield  near  the  fin,  are 
presented.  Experimental  and  numerical  results  are  used  to  describe  the  flowfield  around 
the  single- WAF  configuration  under  consideration.  The  majority  of  the  results  at  Mach  2.8 
have  already  been  published!45’ 108, 109,  no]_ 

First,  the  results  of  pilot  inviscid  calculations  and  photographic  visualization  experi¬ 
ments  are  presented.  Then  results  of  the  viscous  simulation  and  experimental  exploration 
of  the  complex  flowfield  in  the  vicinity  of  the  fin  are  then  presented,  culminating  in  a 
discussion  of  the  aerodynamic  loading. 

5.1  Pilot  Inviscid  Calculations 

The  first  numerical  results  in  this  study  were  provided  by  solving  the  Euler  equations 
over  the  model  geometry  in  the  tunnel’s  test  section  (Figure  5.1).  These  preliminary 
inviscid  results  guided  the  placement  of  the  experimental  measurement  stations. 

5.2  Shadowgraph  and  Schlieren  Photography 

To  get  an  initial  understanding  of  the  flowfield  ahead  of  the  fin,  shadowgraph  and 
schlieren  images  were  obtained  with  the  model  in  the  test  section.  A  shadowgraph  image 
of  the  region  near  the  fin  is  presented  in  Figure  5.2.  Note  how  the  bow  shock  remains 
detached  over  the  full  height  of  the  fin  at  this  Mach  number.  The  compression  waves 
induced  by  the  blended  forebody  of  the  model  can  be  seen  coalescing  into  a  shock  and 
reflecting  off  of  the  opposing  tunnel  wall  just  above  the  fin.  The  bow  shock  created  by  the 
fin  is  also  visible,  although  the  structure  at  the  fin/body  intersection  is  somewhat  obscured 
by  the  boundary  layer  on  the  model  body.  The  boundary  layers  on  the  model  and  opposing 
wall  are  also  identifiable. 
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Figure  5.1  Computed  inviscid  test  section  surface  pressures  (at  M=2.9). 

Composite  schlieren  images  of  the  test  section  with  the  model  installed  are  presented 
in  Figures  5.3  and  5.4.  These  were  made  with  the  horizontal  knife  edge  blocking  the 
top  and  bottom  of  the  light,  respectively,  thus  are  theoretically  negatives  of  one  another. 
The  turbulent  boundary  layer  on  the  model  is  quite  visible  in  Figure  5.3  whereas  the  shock 
structure  near  the  fin/body  juncture  is  most  easily  seen  in  Figure  5.4.  The  blending  shock 
is  again  visible,  as  is  the  bow  shock  created  by  the  fin.  The  other  structures  in  the  images 
which  extend  from  the  floor  upward  are  compression  weak  compression  waves.  These 
disturbances  are  caused  by  small  imperfections  in  the  tunnel  floor  associated  with  the 
removable  plugs  which  are  inserted  into  the  probe  access  slots  and  extend  across  most  of 
the  test  section.  Although  these  disturbances  appear  distinctly,  subsequent  measurements 
have  shown  them  to  be  very  small  compared  to  the  blending  shock  and  the  bow  shock. 
Their  strength  is  exaggerated  in  the  photographs  primarily  due  to  their  two-dimensional 
nature.  In  contrast,  the  bow  and  blending  shocks  are  three-dimensional  structures. 
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Figure  5.2  Shadowgraph  of  fin  region. 


5.3  Lambda  Shock 

Figures  5.5  and  5.6  are  images  produced  by  schlieren  and  shadowgraph  photogra¬ 
phy  which  have  been  enlarged  to  show  the  detail  of  the  Sow  structure  in  front  of  the 
base  of  the  fin.  The  shock-boundary  layer  interaction  produces  the  same  type  of  A-shock 
typically  observed  in  front  of  blunt  fins  and  cylinders  mounted  on  flat  plates  in  super¬ 
sonic  fiowfields.f29’  ^  The  A-shock  structure  was  observed  to  shift  slightly  in  position  from 
photograph  to  photograph. 

Configurations  possessing  blunt  leading  edges  have  a  stagnation  point  on  the  leading 
edge  which  corresponds  to  the  end  of  a  ‘parting  line’  in  the  flowfield,  indicating  the  furthest 
location  from  the  body  at  which  particles  become  entrained  into  the  vortical  flow  ahead  of 
the  leading  edge.  The  investigated  model  has  a  sharp  leading  edge,  making  it  difficult  to 
identify  any  stagnation  point  in  the  numerical  solutions.  However,  particle  traces  of  the 
viscous  numerical  solutions  (Figures,  4.5[c],  5.8[e],&  5.9[e])  indicate  that  the  distance  from 
the  body  at  which  streamlines  diverge  as  they  encounter  the  fin  leading  edge  is  at  yss0.19c. 
This  agrees  quite  well  with  the  A-shock  height  indicated  by  photography  of  y«0.20c. 
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Figure  5.3  Composite  schiieren  photograph  -  knife  edge  on  top. 


Figure  5.4  Composite  schiieren  photograph  -  knife  edge  on  bottom. 
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Figure  5.6  Shadowgraph  of  the  lambda-shock  region. 
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5.4  Effect  of  Fin  Curvature  on  the  Flowfield 

While  2-D  photography  indicates  flow  features  similar  to  those  seen  on  straight  fins, 
conventional  pressure  probes,  hot-film  anemometry,  and  CFD  all  demonstrate  that  the 
flow  near  the  fin  is  highly  asymmetric.!45'  108> 109’  n°l 

It  is  convenient  to  discuss  the  flowfield  in  terms  of  three  regions;  an  upstream  region 
which  is  ahead  of  the  shock  structure,  and  two  downstream  regions,  one  on  either  side  of 
the  fin.  The  present  study  indicates  that  each  of  the  downstream  regions  can  be  further 
divided  into  two  regions;  an  outer  region  which  is  characterized  by  inviscid  behavior,  and 
an  inner  region  near  the  body  in  which  viscous  effects  are  dominant. 

The  flowfield  near  the  single- WAF  geometry  has  been  experimentally  explored!45- 109] 
by  extensively  probing  the  flowfield  near  the  model.  In  addition  to  probing  the  outer  flow 
at  several  locations  near  the  fin  on  the  ceiling  mounted  model!109),  measurements  were 
obtained  nearer  to  the  fin  by  mounting  the  model  on  the  tunnel  side-walls.  The  fuselage 
boundary  layer  was  also  surveyed  at  four  locations  on  the  ceiling  mounted  model  as  shown 
in  Figure  5.7. 145)  These  locations  were  chosen  to  represent  the  upstream  and  downstream 
regions  on  either  side  of  the  fin,  since  two  of  the  stations  set  the  reference  for  the  flow 
upstream  of  the  fin  bow  shock  (at  x=— 0.41c)  and  the  other  two  stations  were  positioned 
downstream  of  the  shock  (at  ®=+0.69c).  At  each  of  these  axial  locations,  the  flow  was 
surveyed  on  the  concave  (Cc)  side  and  on  the  convex  (Cv)  side  of  the  fin  (at  z=±0.47c) 
with  pressure  probes  and  hot-film  cross-wires.  Results  from  the  companion  numerical 
study  were  compared  with  the  experimental  data  at  these  locations,  and  the  combined 
numerical  and  experimental  information  were  examined  for  the  purpose  of  characterizing 
the  flowfield.  Note  that  the  pressure  and  hot-film  probes  flexed  slightly  («2°)  during  the 
experiments  under  aerodynamic  loading.  The  locations  at  which  the  numerical  solution 
was  ‘probed’  were  canted  to  mimic  this  flexing. 

In  the  data  presentation,  the  probe  position  ( x,y,z )  is  nondimensionalized  by  the  fin 
chord,  c=20.3mm,  where  the  coordinate  origin  is  located  at  the  intersection  of  the  body 
surface  centerline  and  leading  edge  of  the  fin.  Negative  x  values  are  upstream  of  the  leading 
edge,  and  negative  z  values  are  to  the  concave  side  of  the  fin.  Boundary  layer  data  are 
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Figure  5.7  Probe  locations. 

presented  as  a  function  of  the  distance,  from  the  model  body,  Y .  This  relative  position 
from  the  body  is  normalized  by  a  reference  boundary  layer  thickness,  do=6-lmm,  which 
was  measured  on  the  model  centerline  0.41c  ahead  of  the  leading  edge  of  the  fin. 

5.4-1  Flow  Ahead  of  the  Bow  Shock.  At  the  upstream  measurement  location,  both 
computed  and  experimental  pitot  pressures  and  mass-flux  profiles  (see  Figures  5.8  and  5.9) 
correspond  to  those  of  a  largely  “undisturbed”  boundary  layer.  While  it  appears  that  the 
numerical  solution  predicts  a  thinner  boundary  layer  than  measured  experimentally,  the 
agreement  is  good  in  the  outer  flow  region.  The  calculations  also  suggest  a  high  degree  of 
flow  symmetry  in  the  outer  flow  at  the  upstream  measurement  locations,  while  the  degree 
of  measured  asymmetry  was  within  the  experimental  uncertainty  (Appendix  C).  At  these 
upstream  locations,  the  numerical  results  indicate  that  the  flow  in  the  boundary  layer  is 
moving  sightly  away  from  the  center!ne(Figures  5.8[d]  and  5.9[d|).  The  flow  very  near 
the  body  is  being  swept  away  from  the  centerline  at  a  very  high  angle,  indicating  this  part 
of  the  boundary  layer  feds  the  presence  of  the  fin. 
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fe]  Calculated  limiting  surface  streamlines,  [f]  Oil  flow  at  Macli=2.06  (Abate  and  Berner) 


Figure  5.9  Numerical  and  experimental  flow  variables;  concave  side  of  fin. 
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[a]  Computed  pressure  contours. 


[b]  Computed  mass-flux  streamlines. 


Figure  5.10  Convex  measurement  plane  (s=-j-0.47c)  pressure  contours  and  mass-flux 
streamlines  ( pu,pv )  given  by  numerical  simulation.  Outline  of  fin  is  over¬ 
laid.  Dashed  lines  represent  survey  locations  in  the  Mach  3  tunnel. 


5.4.2  Flow  on  the  Convex  Side  of  the  Fin.  As  the  flow  nears  the  fin  on  the 
convex  side,  the  outer  flow  (Y/6q>1.5)  passes  through  a  strong  shock  (Figure  5.10).  This 
shock  induces  a  strong  compression  and  deceleration.  As  the  fluid  passes  the  fin  it  is 
expanded  through  a  large  region  of  favorable  pressure  gradient  between  the  shock  and  the 
downstream  measurement  location  (seen  in  Figure  5.10)  due  primarily  to  the  convex  fin 
curvature. 

At  the  downstream  survey  location,  the  pitot  pressure  and  mass-flux  in  the  outer  flow 
have  been  decreased  on  the  convex  side  relative  to  the  upstream  reference  plane.  This  effect 
has  been  captured  by  both  in  viscid-109!  and  viscous  numerical  results  (Figures  5.8[a]&[b]). 
At  this  measurement  station  the  outer  flow  is  directed  away  from  the  body,  but  only 
mildly  away  from  the  fin.  This  behavior  is  seen  in  both  the  experimental  data  and  the 
computational  results  (Figures  5.8[a]&[bj). 

Profiles  of  pressure  and  momentum  in  the  boundary  layer  at  this  location  are  char¬ 
acterized  by  a  large  inflection  (Figure  5.8[a]&[b]).  The  flow  near  Y/S0&  1.1  has  passed  over 
the  horseshoe  vortex  system  produced  by  the  shock/boundary-layer  interaction  ahead  of 
the  fin  (see  Figure  5.10).  In  this  process,  the  flow  greatly  expands  while  only  slightly 
accelerating:  the  net  result  is  a  decrease  in  the  mass-flux.  Flow  in  this  region  is  directed 
strongly  toward  the  body  as  indicated  by  the  inflection  in  the  horizontal  flow  angularity 
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[a]  Numerical  pitot  pressure  and  streamlines.  [b]  Experimental  pitot  pressure. 


Figure  5.11  Flow  at  z=0.69c  measurement  plane  given  by  numerical  simulation  and 
experiment. 


(Figure  5.8[c]).  Agreement  with  experimental  data  is  considered  excellent,  although  the 
tow  turning  angle  is  somewhat  underpredicted.  Examination  of  the  numerical  results  indi¬ 
cate  that  this  turning  effect  is  due  to  a  vortex  embedded  in  the  fin /fuselage  juncture  which 
entrains  fluid,  pulling  it  toward  the  body  (seen  in  Figure  5.11(a)).  This  vortex  creates  the 
pressure  minimum  seen  both  numerically  and  experimentally  in  Figure  5.8(a).  At  roughly 
the  same  location,  a  mild  secondary  flow  component  toward  the  fin  (o<0°)  is  observed  in 
the  numerical  solution;  while  a  corresponding  inflection  is  observed  in  the  experimental 
results  (Figure  5.8(d)).  Although  the  numerical  and  experimental  6  profiles  are  somewhat 
different,  the  discrepancies  may  be  attributed  to  probe  volume  effects  and  probe  location 
errors  (see  Section  C.4).  Notably,  if  the  location  at  which  the  numerical  solution  is  ex¬ 
amined  is  shifted  by  the  probe  volume  (1mm)  closer  to  the  fin,  the  predicted  inflection 
in  4>  very  closely  resembles  the  experimentally  obtained  data  (Figure  5.12(b)).  Also,  the 
flattening  of  the  B  profile  in  the  experimental  data  over  the  range  0.3<Y/£o<0.8  is  closely 
duplicated  (Figure  5.12(a)).  It  is  notable  that  the  juncture  vortex  and  the  associated 
acceleration  toward  the  body  are  viscous  phenomena  and  thus  have  not  been  captured  by 
inviscid  methods.!109!  The  significance  of  this  finding  is  addressed  in  Section  5.5. 

Slightly  closer  to  the  body  (F/<5oc[0.3,1.0]),  the  flow  experiences  a  compression  from 
above  while  is  at  the  same  time  aligned  with  the  z-axis  near  the  body.  The  net  effect  is 
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[a]  Horizontal. flow  angularity,  6  (degrees).  [b]  Azimuthal  flow  angularity,  <f>  (degrees). 

Figure  5.12  Effect  of  Az=-lmm  on  flow  angularity;  convex  side  of  fin. 

a  sharp  increase  in  mass-flux.  Below  F/^o«0.3,  wall  effects  force  a  decrease  in  mass-flux 
and  pitot  pressure.  The  flow  is  directed  downward  and  away  from  the  fin  over  a  very  small 
region  (F/<5o<0.2),  following  the  contour  of  the  body.  The  numerical  results  and  oil  flow 
patterns  at  Mach  2.06 I2)  (discussed  below)  suggest  that  the  azimuthal  flow  angularity,  <j>, 
at  this  location  tends  toward  zero  at  the  wall. 

Given  its  proximity  to  the  fin,  and  hence  its  effects  on  the  aerodynamic  loading, 
more  discussion  on  the  ‘juncture’  vortex  is  warranted.  The  juncture  vortex  (Figure  5.13) 
originates  near  the  leading  edge  of  the  fin/body  juncture  and  remains  tucked  into  the 
fin/body  junction,  growing  in  strength  and  size  as  it  progresses  along  the  fin.  The  size 
and  orientation  of  this  vortex  is  clearly  evident  in  limiting  surface  streamlines  calculated 
from  the  numerical  solution  (Figure  5.8[e])  and  in  the  surface  oil  flow  patterns  obtained 
by  Abate  and  Berner!2]  at  Mach  2.06  (Figure  5.8[f]).  Surface  streamlines  starting  at  the 
leading  edge  travel  downward  along  the  beveled  edge  and  join  with  streamlines  flowing  up 
from  the  root  to  form  an  accumulation  of  oil  film  (or  convergence  of  streamlines)  on  the 
surface.  This  convergence  line  marks  the  separation  line  formed  by  the  juncture  vortex, 
and  moves  away  from  the  juncture  as  it  travels  toward  the  trailing  edge.  The  complicated 
flow  structure  observed  in  the  oil  flow  patterns  closely  resembles  that  predicted  by  the 
numerical  solution,  suggesting  that  the  flow  structure  near  the  juncture  changes  little 
within  this  Mach  number  range.  On  the  body,  a  weak  attachment  line  (surface  streamline 
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Figure  5.13  Stagnation  pressure  iso-surfaces  and  streamlines  in  juncture  vortex  region  on 
convex  side  of  the  fin. 

divergence)  moving  outward  from  the  leading  edge  is  clearly  evident  in  both  the  numerical 
solution  (Figure  5.8[ej)  and  the  oil  flow  (Figure  5.8[f]). 

Evidence  of  such  a  vortex  has  also  been  observed  in  oil  flow  patterns  on  straight  blunt 
fins  mounted  on  flat  plates  J24-  Such  vortices  have  been  observed  to  change  rotational 
direction  on  straight  fins  depending  on  incidence  angle.  The  viscous  numerical  results 
indicate  that  the  rotation  of  the  juncture  is  of  the  same  sense  as  that  seen  on  the  com¬ 
pression  side  of  a  straight  fin  at  incidence.  Thus,  with  respect  to  the  juncture  vortex,  fin 
curvature  and  attachment  angle  can  induce  similar  effects  to  those  produced  by  cross-flow. 


5-4-3  Flow  on  the  Concave  Side  of  the  Fin.  In  contrast  to  the  flow  on  the  convex 
side  of  the  fin,  the  flow  on  the  concave  side  passes  through  a  somewhat  weaker  shock 
(Figure  5.14).  Thus,  the  flow  undergoes  a  much  more  modest  deceleration.  Also,  the 
post-shock  expansion  is  partially  offset  by  the  compressive  effects  of  fin  curvature.  The  net 
effect  is  a  dramatic  increase  in  the  mass-flux  (up  to  30%)  at  the  downstream  measurement 
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[a]  Computed  pressure  contours. 


[b]  Computed  mass-flux  streamlines. 


Figure  5.14  Concave  measurement  plane  (z= -0.47c)  pressure  contours  and  mass-flux 
streamlines  (pu,pv)  given  by  numerical  simulation.  Outline  of  fin  is  overlaid. 
Dashed  lines  represent  survey  locations  in  the  Mach  3  tunnel. 


location  as  compared  to  the  upstream  location;  This  is  evident  in  both  the  numerical  and 
experimental  data.  On  this  side  of  the  fin,  the  outer  flow  is  strongly  directed  away  from 
the  fin  (Figure  5.9[d])  at  flow-  angles,  O.  up  to  10®  at  the  mid-span  (y/<!>o~2.5).  Here, 
numerical  and  experimental  results  indicate  that  0=0°,  meaning  that  the  flow  is  directed 
toward  the  center  of  fin  curvature  (Figures  5.9[c]  and  5. Ufa]). 

Approaching  the  body,  the  fluid  momentum  decreases  (Figure  5.9[al&[b]).  Over  a 
small  region  inside  the  boundary  layer  (F/^fO.o,  1.1]),  the  numerical  solution  suggests 
that  there  is  a  large  inflection  in  the  azimuthal  angularity  (Figure  5.9[d])  where  the  flow  is 
almost  aligned  with  the  vertical  plane  (<£» 0°).  This  inflection  is  more  pronounced  in  the 
experiment  data,  but  occurs  at  the  same  location.  This  effect  is  likely  to  be  a  combination 
of  the  flow  wrapping  around  the  fin  and  an  expansion  which  reflects  off  of  the  bow  shock  as 
a  compression.  Inviscid  numerical  results!45'  109, 1103  on]y  faintly  hinted  at  this  trend.  Over 
this  same  range,  the  experimental  data  (Figure  5.9[c])  suggest  that  the  magnitude  of  the 
horizontal  flow  angularity  is  greatly  reduced  (i.e.  6—0  at  Y/$o~0.Z).  However,  this  change 
in  flow  angularity  was  not  present  in  the  numerical  solution  at  the  nominal  measurement 
location.  This  discrepancy  can  be  attributed  to  probe  volume  effects  and  probe  location 
errors  (see  Section  C.4).  Notably,  this  trend  is  clearly  present  in  the  numerical  solution  at 
a  probe  width  (1mm)  further  from  the  fin  (Figure  5.15[a]). 
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[a]  Horizontal  flow  angularity,  9  (degrees).  [b]  Azimuthal  flow  angularity,  <j>  (degrees). 

Figure  5.15  Effect  of  Az=+lmm  on  flow  angularity;  concave  side  of  fin. 

Closer  to  the  body  (Y/%<0.5)  there  is  a  small  region  in  which  measured  and  com¬ 
puted  pitot  pressures  do  not  change.  The  numerical  results  suggest  that  the  flow  is  moving 
downward  and  away  from  the  fin  (0«-3O°,<£«— 39°),  lending  to  the  flattening  of  pitot 
pressure.  Cross- wire  volume  effects  precluded  experimental  examination  of  this  region.  (451 

On  this  side  of  the  fin,  surface  streamlines  starting  at  the  leading  edge  also  travel 
downward  along  the  beveled  edge  and  join  with  streamlines  flowing  up  from  the  root. 
However  the  streamline  convergence  is  ‘incomplete’  from  below,  and  no  juncture  vortex  is 
indicated.  As  on  the  convex  side,  the  similarities  between  the  predicted  surface  streamlines 
and  the  observed  oil  flow  patterns  of  Abate  and  Berner  at  Mach  2.06 12!  suggest  that  the 
flow  structure  near  the  juncture  on  this  side  of  the  fin  changes  little  within  this  Mach 
number  range.  The  downstream  measurement  station  on  this  side  of  the  fin  is  located 
just  behind  a  separation  line  on  the  body,  which  is  seen  in  the  computational  results  and 
oil  flow  pattern.  Due  to  the  oblique  attachment  angle  («135°),  no  juncture  vortex  was 
observed  on  the  concave  side  of  the  fin  in  either  the  numerical  or  experimental  studies,  nor 
is  one  indicated  by  the  oil  flow  visualizations  I2 1  at  Mach  2.06  (Figures  5.9[e]&[f]). 

In  comparing  the  numerical  results  to  the  oil  flow  patterns  of  Abate  and  Berner  t2!, 
it  is  notable  that  the  latter  was  obtained  on  a  four-finned  missile.  Thus  the  similarity 
between  fin  surface  streamline  patterns  suggests  that  the  single-fin  model  produces  the 
relevant  flow  features  present  on  configurations  with  multiple  wrap-around  fins. 
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5.5  Aerodynamic  Loading  on  the  Fin 

The  diverse  flow  topologies  on  either  side  of  the  WAF  produce  an  asymmetric  load 
distribution  on  the  fin  surface.  While  previous  inviscid  calculations  have  captured  many  of 
the  essential  flow  features,  the  juncture  vortex  on  the  convex  side  of  the  fin  is  dominated 
by  viscous  effects  and  its  position  and  growth  near  the  fin  root  may  provide  a  significant 
aerodynamic  load,  particularly  at  higher  Mach  numbers. 

The  fin  surface  pressures  computed  in  the  inviscid  and  viscous  simulations  are  shown 
in  Figures  5.16  and  5.17,  respectively.  On  the  concave  side  of  the  fin,  the  bow  shock  is 
“focused”,  producing  a  high  pressure  region  between  the  fin  and  its  center  of  curvature 
(Figure  5.11)  where  mass  flux  levels  are  increased  to  50%  over  the  free-stream  value.  This 
produces  large  region  of  relatively  high  surface  pressures  near  the  half-span  of  the  fin 
that  contributes  to  a  negative  rolling  moment  (Figures  5.16[b]  and  5.17[b]).  Here,  rolling 
moment  will  be  defined  in  the  vehicle  stability  sense,  thus  a  negative  value  indicates  a 
moment  acting  in  the  direction  of  negative  curvature.  The  inviscid  calculations  predicted 
that  the  convex  side  of  the  fin  also  had  a  region  near  the  fin  root  over  which  high  pressure 
levels  nearly  reach  the  magnitudes  seen  on  the  concave  side.  This  compression  had  been 
attributed  to  the  fin  being  canted  in  the  convex  direction  (w45°)  at  the  fin/body  inter¬ 
section.!109)  However,  the  viscous  numerical  results  show  that  this  high  pressure  region  is 
displaced  by  the  boundary  layer,  and  weakened  greatly  (Figure  5.17[a]).  Instead  of  high 
pressure,  the  root  region  is  dominated  by  low  pressures  induced  by  the  juncture  vortex. 
Using  inviscid  numerical  methods,  Abate  and  Cook!3)  have  shown  that  the  rolling  moment 
is  a  function  of  both  the  fin  curvature  and  fin  attachment  angle.  However,  it  is  now  clear 
that  the  effect  of  fin  attachment  angle  is  not  fully  captured  by  an  inviscid  analysis.  While 
the  Euler  analysis  predicted  a  rolling  moment  coefficient  of  —0.0102,  the  viscous  simulation 
predicted  a  value  of  -0.0112,  a  10%  increase. 
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5.6  Effect  of  the  Fin  on  Fuselage  Boundary  Layer  Turbulence 

Turbulence  quantities  were  also  measured  at  the  same  locations  discussed  in  Sec¬ 
tion  5.4  (see  Figure  5.7).  t4^!  Although  the  downstream  probe  location  on  the  concave  side 
of  the  fin  is  in  a  region  of  mild  favorable  pressure  gradient,  it  is  also  very  close  to  the 
bow  shock.  Examination  of  the  turbulence  intensity  in  this  region  indicates  that  the  shock 
induces  an  increase  in  the  axial  (Figure  5.18[a])  and  transverse  (Figures  5.18[b]  and  5.18[c]) 
turbulence  intensity,  which  is  consistent  with  the  destabilizing  effects  of  adverse  pressure 
gradients  on  turbulent  boundary  layers. ("1 

On  the  convex  side,  the  turbulence  intensities  at  the  same  streamwise  location  are 
far  lower,  since  the  flow  has  experienced  a  favorable  pressure  gradient  over  an  extended 
streamwise  distance  (Figure  5.10).  Expansions  have  been  reported  to  stabilize ,  or  reduce, 
the  turbulence  levels.  Also  contributing  to  this  dramatic  recovery  is  the  rapid  flow 
acceleration  induced  by  the  convex  curvature  of  the  fin.  Both  of  the  secondary  mass- 
flux  turbulence  intensities  experience  a  sharp  rise  near  Y/6 o=l.l,  but  not  so  the  axial 
component.  The  rise  in  cross-flow  turbulence  intensity  occurs  at  the  same  location  where 
the  mean  flow  is  being  turned  sharply  toward  the  body  by  the  juncture  vortex. 

The  net  effect  of  the  bow  shock  and  fin  curvature  are  illustrated  in  Figures  5.19 
and  5.20  via  the  nondimensional  turbulent  kinetic  energy  (TKE),  Kc,  defined  by 


The  TKE  is  significantly  elevated  on  the  concave  side  of  the  fin  and  reduced  on  the  convex 
side  relative  to  values  upstream  of  the  shock  (Figure  5.19). 

Figure  5.20  shows  the  nondimensional  TKE  at  two  streamwise  stations;  one  at  the 
same  location  as  the  boundary  layer  measurements  (x=+ 0.69c),  and  one  further  upstream 
at  x=+0.38c.  These  surveys,  which  included  the  upper  portion  of  the  boundary  layer, 
indicate  that  TKE  levels  on  the  convex  side  are  markedly  greater  than  those  on  the  concave 
side  at  both  stations.  Also,  TKE  is  dissipating  in  the  axial  direction  on  both  sides  of  the 
fin,  as  the  flow  passes  through  regions  of  favorable  pressure  gradient.  Thus,  as  the  flow 
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Figure  5.18  Turbulence  intensity  profiles 
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Figure  5.20  Turbulent  kinetic  energy,  Kc. 


continues  to  recover,  the  TKE  levels  on  the  concave  side  of  the  fin  are  likely  to  continue 
decreasing  in  the  downstream  direction,  possibly  to  the  levels  seen  on  the  convex  side. 

The  stabilizing  and  destabilizing  effects  on  turbulence  to  either  side  of  the  fin  are 
also  revealed  in  the  Reynolds  shear  stress  estimates  measured  with  the  cross- wire  (Equa¬ 
tion  3.6,  Chapter  3)  presented  in  Figure  5.21.  Upstream  of  the  fin,  the  shear  stress  profile 
corresponds  to  that  of  an  undisturbed  boundary  layer,  with  levels  comparable  to  those 
upstream  of  the  model.  This  indicates  that  the  effects  of  the  compression  caused  by  the 
blended  region  of  the  model  have  been  damped  to  levels  comparable  to  an  equilibrium 
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Figure  5.21  Measured  turbulent  shear  stresses,  r?  =  _|_  p  Ui  u.  ^40  . 

turbulent  boundary  layer.  Note  that  the  second  term  of  the  total  turbulent  shear  stress  in 
Equation  3.6  has  been  determined  to  be  at  least  an  order  of  magnitude  smaller  than  the 
first  term  for  all  surveyed  regions,  even  in  those  having  large  flow  angularity.  Therefore, 
only  the  total  turbulent  shear  stress  is  presented. 

On  the  convex  side  of  the  fin,  a  reduction  of  turbulent  shear  stress  with  downstream 
position  is  indicative  of  the  strong  stabilizing  effect  of  favorable  pressure  gradient.  Indeed, 
as  observed  by  other  researchers  investigating  the  correlation  between  streamline  distortion 
and  turbulence!15!,  the  expansion  associated  with  a  favorable  pressure  gradient  can  result 
in  reduced  (or  even  negative)  turbulent  shear  stress.  Conversely,  the  shear  stress  on  the 
concave  side  is  increased  by  roughly  100-200%,  commensurate  with  the  previously  noted 
increases  in  turbulence  intensity  and  turbulent  kinetic  energy.  This  pattern  is  consistent 
with  other  measurements  obtained  in  regions  of  large  compression,  t15!  Note  too  that  the 
boundary  layer  thickness  has  dramatically  increased  (by  about  60-70%)  on  the  concave 
side,  while  it  was  reduced  slightly  on  the  convex  side.  The  large  increase  in  turbulent 
shear  stress  on  the  concave  side  of  the  fin  may  be  the  primary  reason  that  measured  pitot 
pressure  and  mass  flux  levels  were  less  than  those  predicted  by  the  numerical  simulation  in 
the  outer  boundary  layer  (Figure  5.9[a]&[b]),  since  the  Baldwin-Lomax  turbulence  model 
was  not  designed  to  simulate  the  effects  of  pressure  gradient  on  turbulence. 
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6.  WAF  Flowfield  at  Mach  4-9 

The  single- WAF  model  has  also  been  investigated  experimentally  and  numerically  at 
Mach  4.9  and  a  unit  Reynolds  number  of  Re/^«65-75xl06m-1.  These  investigations  were 
aimed  at  determining  the  structure  of  the  mean  flowfield  at  a  high-speed  condition. 

In  the  experiment,  the  mean  flow  around  the  ceiling-mounted  WAF  configuration 
was  surveyed  using  pitot  and  cone-static  probes  and  was  visualized  using  shadowgraph  and 
schlieren  photography  in  the  AFIT  Mach  5  wind  tunnel  (Appendix  D).  Fuselage  boundary 
layer  surveys  were  obtained  upstream  and  downstream  of  the  bow  shock.  The  fuselage 
boundary  layer  flow  was  explored  near  the  fin  on  the  ceiling  mounted  model,  as  it  was  in  the 
Mach  3  tunnel.  Placement  of  the  measurement  stations  for  the  experiment  was  guided  by 
the  inviscid  computations  (Figure  6.1)  and  by  experience  gained  from  the  experiment  in  the 
Mach  3  tunnel.  The  probe  locations  were  the  same  in  the  span- wise  direction  (z=±0.47c) 
as  they  were  in  the  Mach  2.8  experiment,  and  the  stations  that  set  the  reference  for  the 
flow  were  again  placed  just  upstream  of  the  flow  interaction  at  x= -0.41c.  However,  the 
probes  were  positioned  further  aft  at  the  downstream  locations  (a;=+0.84c)  so  as  to  remain 
well  behind  the  bow  shock. 

In  the  viscous  numerical  simulation,  the  Navier-Stokes  equations  were  solved  using 
the  same  computational  grid  and  numerical  strategy  used  for  the  Mach  2.8  simulation.  The 
algebraic  turbulence  model  of  Baldwin  and  Lomax  I11!  was  again  employed.  Experimental 
pressure  data  were  used  to  define  the  upstream  boundary  condition  for  this  calculation. 
The  resulting  numerical  solution  is  compared  with  the  experimental  data  and  the  combined 
sets  of  information  are  examined  to  characterize  the  flowfield. 

In  the  presentation  of  the  results,  the  probe  position  ( x,y,z )  is  nondimensionalized 
by  the  fin  chord  as  it  was  for  the  Mach  2.8  results.  Boundary  layer  data  are  presented  as 
a  function  of  the  distance  from  the  model  body,  Y.  This  relative  position  from  the  body 
is  normalized  by  a  reference  boundary  layer  thickness,  <5co=10.2mm,  which  was  measured 
on  the  tunnel  centerline  0.8r  ahead  of  the  blended  body.  Note  that  this  is  larger  than 
the  dimension  (<5o=6.2mm)  by  which  the  Mach  2.8  results  were  normalized.  Although 
the  pressure  probes  used  in  the  Mach  5  wind  tunnel  experienced  higher  loadings  than 
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Figure  6.1  Computed  inviscid  test  section  surface  pressures  (at  M=5.0). 

the  probes  used  in  the  Mach  3  tunnel,  they  were  much  stiffer  and  experienced  less  flexing 
(«1.1°).  The  locations  at  which  the  numerical  solution  is  compared  to  experimental  results 
are  canted  to  mimic  this  flexing. 


6.1  Shadowgraph  and  Scklieren  Photography 

Shadowgraphs  and  schlieren  images  of  the  fin  region  (Figures  6.2  and  6.3)  again 
indicate  that  fin  shock  remained  detached  over  the  full  span  of  the  fin  as  it  did  at  Mach  2.8, 
although  stand-off  distance  was  reduced.  The  same  principal  features  observed  at  the 
lower  Mach  number  were  again  visible,  including  the  shock  caused  by  the  blended  forebody 
and  the  bow  shock,  both  of  which  v/ere  more  highly  swept  than  at  Mach  2.8.  The  A-shock 
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Figure  6.3  Schlieren  photograph  (M=4.9). 


6-3 


[a]  x=0.69c. 


[b]  x=0.B4c. 


Figure  6.4  Computed  pitot  pressure  and  secondary  streamlines  at  M=2.8  measurement 
plane,  and  at  present  measurement  plane  (a:=0.84c). 


was  also  distinct,  though  somewhat  unsteady,  and  positioned  slightly  closer  to  the  body 
and  at  a  shallower  angle  than  observed  in  the  Mach  2.8  experiment  (c/ Figures  5.2  and  6.2). 
This  trend  was  also  captured  in  the  viscous  numerical  solutions  (c/Figures  4.5  and  4.8). 


6.2  Effect  of  Fin  Curvature  on  the  Mean  Flowfield 

The  experimental  and  numerical  results  suggest  that  the  flow  near  the  fin  is  highly 
asymmetric  as  it  was  at  Mach  2.8.  In  fact,  most  of  the  qualitative  discussion  on  the 
flowfield  structure  at  Mach  2.8  (Chapter  5)  applies  at  Mach  4.9  as  well.  As  expected,  the 
fin’s  domain  of  influence  in  the  outer  flow  w-as  reduced  from  that  at  Mach  2.8  (Figure  6.4). 


The  measured  flow  asymmetry  at  the  upstream  measurement  plane  was  minimal 
(Figure  6.5),  and  well  within  the  experimental  uncertainty  (especially  considering  that  the 
data  were  acquired  on  different  days).  At  this  upstream  location,  there  is  a  much  larger 
pressure  gradient  in  the  direction  away  from  the  body  than  was  observed  at  Mach  2.8  due 
to  the  relative  proximity  of  the  blending  shock. 
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Figure  6.5  Experimental  and  numerical  pitot  pressure. 

6.2.1  Flow  on  the  Convex  Side  of  the  Fin.  As  seen  at  Mach  2.8,  the  flow  on 
the  convex  side  of  the  fin  passes  through  a  strong  shock  (Figure  6.6[a]),  and  then  expands 
through  a  large  region  of  favorable  pressure  gradient  between  the  shock  and  the  downstream 
measurement  location  due  to  the  convex  fin  curvature.  Relative  to  its  upstream  value, 
the  pitot  pressure  in  the  outer  flow  at  the  downstream  measurement  station  has  been 
decreased  by  about  40%  (Figure  6.5[a]).  Agreement  between  experimental  and  computed 
pitot  pressures  is  considered  excellent  on  this  side  of  the  fin. 

The  calculated  pitot  pressure  profile  features  a  large  inflection  in  the  inner  region 
(Figure  6.5[a]).  This  inflection  is  more  difficult  to  identify  in  the  experimental  results  than 
it  was  at  Mach  2.8.  While  the  larger  volume  of  the  pitot  probe  used  at  Mach  5  precluded  a 
definitive  assessment  of  this  feature,  it  appears  that  the  inflection  may  be  closer  to  the  body 
and  not  as  large  as  that  predicted  by  the  numerical  simulation.  The  computed  pressure 
peak  is  at  approximately  the  same  physical  distance  from  the  wall  at  both  Mach  numbers 
(0.3*00*0.5*0). 

Recall  from  the  Mach  2.8  discussion  that  this  inflection  in  pitot  pressure  was  asso¬ 
ciated  with  the  existence  of  a  juncture  vortex.  This  vortex  is  also  present  in  the  solution 
at  Mach  4.9  (Figure  6.7),  and  is  probably  present  in  the  experimental  flowfield,  although 
the  pressure  data  do  not  provide  conclusive  evidence.  The  predicted  rotational  direction 
at  Mach  4.9  is  the  same  as  it  was  at  Mach  2.8,  but  the  predicted  location  of  the  junc- 
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Figure  6.6  Computed  pressure  and  mass-flux  streamlines  (pu,pv)  on  measurement  planes 
(M=4.9).  Outline  of  fln  is  overlaid.  Dashed  lines  represent  survey  locations. 


Figure  6.7  Stagnation  pressure  iso-surfaces  and  streamlines  in  juncture  vortex  region  on 
convex  side  of  the  fin. 
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Figure  6.8  Computed  limiting  surface  streamlines  at  M=4.9. 

ture  vortex  is  slightly  closer  to  the  fuselage  (cf  Figures  5.11  [a]  and  6.4[a]).  The  size  and 
orientation  of  this  vortex  is  also  indicated  by  the  calculated  limiting  surface  streamlines 
(Figure  6.8fa]).  While  the  structure  of  the  surface  streamlines  on  the  fin  are  very  similar  to 
those  calculated  at  Mach  2.8,  the  separation  line  formed  by  the  juncture  vortex  is  slightly 
closer  to  the  body.  Also,  the  weak  attachment  line  on  the  body  moving  outward  from  the 
leading  edge  is  slightly  closer  to  the  fin  than  predicted  at  Mach  2.8. 

6.2.2  Flow  on  the  Concave  Side  of  the  Fin.  As  was  the  case  at  Mach  2.8,  the 
outer  flow  (F/doo>1.0)  on  the  concave-side  measurement  plane  passes  through  a  somewhat 
weaker  shock  than  on  the  convex  side  (Figure  6.6[bj).  Again,  the  post-shock  expansion  is 
partially  offset  by  the  compressive  effects  of  fin  curvature.  The  net  effect  is  a  very  large 
increase  in  the  pitot  pressure  (up  to  90%)  and  momentum  by  the  time  the  flow  reaches 
the  downstream  measurement  location.  This  increase  is  observed  in  both  the  numerical 
and  experimental  data  (Figure  6.5[b]),  and  is  even  more  dramatic  than  at  Mach  2.8.  The 
overprediction  of  the  pitot  pressure  in  the  outer  boundary  layer  on  this  side  of  the  fin  is 
probably  again  the  result  of  the  lack  of  pressure  gradient  effects  in  the  turbulence  model. 

Close  to  the  body  (Y/6oo<0A)  there  is  a  small  region  over  which  measured  and 
numerical  pitot  pressures  do  not  change.  The  viscous  simulations  suggest  that  the  flow  in 
this  region  is  moving  down  and  away  from  the  fin  at  angles  comparable  to  those  predicted 
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Figure  6.9  Computed  inviscid  fin  surface  pressures. 


at  Mach  2.8,  leading  again  to  a  ‘flattening’  of  pitot  pressure.  The  similarities  between 
the  calculated  flowfields  the  two  Mach  numbers  (M=2.8,4.9),  as  well  as  the  similarity  of 
calculated  surface  streamline  patterns  to  the  oil  flow  patterns  of  Abate  and  Berner  at 
Mach  2.06®,  suggest  that  the  flow  structure  on  the  concave  side  of  the  fin  changes  little 
■with  Mach  number.  The  exception  is  the  trajectory  of  the  incomplete  surface  streamline 
convergence  on  the  fin,  which  becomes  more  obtuse  "with  decreasing  Mach  number.  The 
downstream  measurement  station  on  this  side  of  the  fin  is  again  located  just  behind  a 
separation  line  on  the  body  which  is  at  approximately  the  same  location  as  it  was  at 
Mach  2.8.  On  this  side  of  the  fin,  the  calculated  secondary  flow  structure  is  also  very 
similar  to  that  predicted  at  Mach  2.8  (c/  Figure  6.4[a]  and  5.11[a]).  The  outer  flow  is 
strongly  directed  away  from  the  fin  toward  the  center  of  fin  curvature  at  the  mid-span, 
and  again,  no  juncture  vortex  is  observed  on  this  side  of  the  fin. 


6.3  Aerodynamic  Loading  on  the  Fin 

As  was  true  at  Mach  2.8,  diverse  flow  topology'  on  either  side  of  the  WAF  produce  a 
dramatically  different  load  distributions  on  the  opposing  fin  surfaces.  The  fin  surface  pres¬ 
sures  predicted  by  the  inviscid  and  viscous  simulations  are  shown  in  Figures  6.9  and  6.1Q, 
respectively.  The  computed  rolling  moment  from  this  simulation  at  Mach  4.9  is  about 
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Figure  6.10  Computed  viscous  turbulent  (Baldwin-Lomax)  fin  surface  pressures. 


one-third  of  that  predicted  at  Mach  2.8,  and  is  about  60%  of  that  predicted  by  inviscid 
simulation.  As  shown  in  Figure  6.11,  the  viscous  simulations  suggest  that  Mach  number 
has  a  stronger  influence  on  rolling  moment  that  inviscid  simulations  would  imply. 


Figure  6.11  Computed  rolling  moments. 
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7.  Conclusions  and  Recommendations 


7.1  Conclusions 

The  structure  of  the  fiowfield  near  a  single  wrap-around  fin  (WAF)  mounted  on 
a  semi-cylindrical  body  has  been  characterized  using  both  experimental  and  numerical 
methods  at  Mach  numbers  of  2.8  and  4.9.  A  single-finned  model  was  first  extensively 
tested  in  the  AFIT  Mach  3  wind  tunnel.  In  this  experiment,  the  flow  around  the  test 
article  was  surveyed  at  several  stations  along  its  length,  concentrating  on  the  region  near 
the  fin.  While  the  boundary  layer  on  the  fin  was  determined  to  be  too  thin  to  survey 
with  pressure  or  hot-film  cross-wire  probes,  the  boundary  layer  on  the  missile  fuselage  was 
easily  explored  using  these  devices.  The  result  was  a  mapping  of  the  pressure,  velocity, 
and  turbulent  properties  near  the  fin.  The  mean  flow  near  the  model  was  also  obtained  in 
the  AFIT  Mach  5  wind  tunnel.  Taken  together,  these  experimental  studies  comprise  a  set 
of  mean  flow  and  turbulence  data  not  previously  available  for  curved  fins. 

Companion  numerical  studies  were  also  performed  wherein  the  Reynolds  averaged 
Navier-Stokes  equations  were  solved  with  the  algebraic  turbulence  model  of  Baldwin  and 
Lomaxl11!  in  the  vicinity  of  the  single- WAF  model.  The  excellent  agreement  with  experi¬ 
mental  data  suggests  that  the  calculations  have  captured  the  relevant  flow  physics  involved 
in  this  complicated  fiowfield.  It  is  notable  that  the  oil  flow  pattern  of  Abate  and  Bernert2! 
to  which  these  results  compared  so  favorably  was  obtained  on  a  four-finned  missile.  Thus 
the  resemblance  of  computed  and  observed  surface  streamline  patterns  on  the  fin  suggests 
that  the  present  study,  which  uses  a  simplified  single-fin  model,  captures  the  relevant  flow 
features  in  the  fin  region  for  a  non-spinning  missile  with  multiple  wrap-around  fins.  Taken 
in  concert,  the  experimental  and  numerical  results  have  been  interpreted  to  characterize 
the  fiowfield  in  the  vicinity  of  a  wrap-around  fin.  Based  on  the  results  of  this  research, 
several  conclusions  may  be  made  regarding  the  nature  of  the  fiowfield  near  wrap-around 
fins. 

One  of  the  more  significant  findings  of  the  present  study  is  that  both  inviscid  and 
viscous  properties  play  significant  roles  in  determining  the  structure  of  the  fiowfield  near 
WAFs.  The  outer  fiowfield  exhibits  asymmetries  brought  about  by  the  effects  of  pressure 
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gradient,  streamline  curvature,  and  differing  shock/expansion  structures  -  while  viscous 
phenomena  induce  asymmetries  near  the  fuselage.  Regarding  the  latter,  the  Navier- Stokes 
simulations  predicted  a  vortex  in  the  fin/body  juncture  on  the  convex  side  of  the  fin.  This 
vortex,  not  present  on  the  concave  side  presumably  due  to  the  oblique  fin  attachment  angle, 
increases  the  pressure  loading  near  the  fin  root.  The  net  result  is  a  pressure  differential 
across  the  fin  which  alters  the  rolling  moment.  Both  hot-film  anemometry  and  surface  flow 
visualizations  corroborate  the  existence  of  this  viscous-induced  vortical  structure.  Inviscid 
simulations  cannot  capture  this  vortex  and  thus  may  not  be  expected  to  reasonably  predict 
the  stability  behavior  of  missiles  having  WAFs. 

That  said,  many  aspects  of  the  flowfield  were  accurately  captured  by  solving  the 
Euler  equations.  The  bow  shock  created  by  the  wrap-around  fin  is  an  inviscid  phenomenon. 
Except  in  the  immediate  vicinity  of  the  fin/body  intersection  where  the  bow  shock  interacts 
with  the  missile  body  boundary  layer,  the  pressure  field  is  dominated  by  inviscid  effects. 
This  has  been  demonstrated  through  the  excellent  agreement  between  measured  quantities 
and  those  predicted  by  inviscid  (as  well  as  viscous)  numerical  methods  in  the  outer  region 
of  the  flow.  The  shock  remains  detached  over  the  full  span  of  the  fin  at  Mach  2.8  and  at 
Mach  4.9,  and  its  interaction  with  the  fuselage  boundary  layer  creates  the  same  type  of 
A-shock  associated  with  blunt  fins  in  supersonic  flowfields. 

The  reduction  of  data  from  the  pitot,  cone-static  and  hot-film  probes  at  Mach  2.8  pro¬ 
duced  a  significant  amount  of  turbulence  data.  Prior  to  this  study,  no  detailed  mean  flow  or 
turbulence  measurements  existed  for  WAF  missile  configurations.  These  data  yielded  some 
interesting  insights.  As  expected,  the  bow  shock  causes  a  dramatic  increase  in  turbulent 
kinetic  energy  and  Reynolds  shear  stress  on  both  sides  of  the  fin.  The  flow  experiences  an 
expansion  as  it  passes  the  fin  which  reduces  the  turbulence  intensity.  However,  for  a  fixed 
streamwise  location,  the  reduction  is  far  greater  on  the  convex  side  of  the  fin,  where  the 
flow  experiences  a  stronger  favorable  pressure  gradient  over  a  longer  distance.  This  results 
in  lower  turbulence  intensities,  producing  lower,  though  still  significant,  shear  stresses.  It 
is  notable  that  the  turbulence  model  used  was  not  designed  to  account  for  the  effects  of 
pressure  gradient  and  streamline  curvature.  This  may  be  largely  responsible  for  the  over- 
prediction  of  momentum  levels  in  the  outer  boundary  layer  on  the  concave  side  of  the  fin  in 


7-2 


the  present  numerical  simulations.  It  is  expected  that  these  turbulence  data  will  be  useful 
for  validation  of  turbulence  closure  models  intended  to  predict  flows  having  large  pressure 
gradients. 

At  Mach  4.9,  both  the  numerical  and  experimental  results  qualitatively  resemble 
those  obtained  at  Mach  2.8.  Examination  of  the  computed  surface  streamline  patterns 
suggests  that  the  flow  structure  in  the  inner  viscous  region  is  somewhat  invariant  over  this 
range  of  Mach  number.  However,  the  juncture  vortex  is  observed  to  shift  slightly  toward 
the  body  with  increasing  Mach  number,  causing  a  dramatic  reduction  in  computed  rolling 
moment  compared  to  that  predicted  by  inviscid  theory.  While  the  viscous  computations 
did  not  predict  a  rolling  moment  reversal  at  either  of  the  Mach  numbers  investigated,  the 
predicted  rolling  moment  had  a  much  greater  dependence  on  Mach  number  (in  a  manner 
consistent  with  a  moment  reversal)  than  inviscid  simulations  have  projected.  Again,  this 
underscores  the  importance  of  viscous  effects  on  the  rolling  moment. 

7.2  Recommendations 

Although  the  present  study  has  provided  clearer  picture  of  the  flowfleld  near  a  WAF 
in  the  given  environment,  it  was  by  necessity  limited  in  scope,  and  thus  represents  a  first 
step  toward  understanding  the  flowfleld  dynamics  of  deployed  missiles  employing  WAFs. 
What  follows  are  suggested  research  areas  that  could  expand  the  understanding  of  WAF 
aerodynamics.  Unquestionably,  this  problem  should  continue  to  be  addressed  using  both 
experimental  and  numerical  means. 

7.2.1  General.  As  a  means  of  addressing  the  problems  induced  by  the  variation  in 
rolling  moment  with  Mach  number,  the  current  research  suggests  that  remedies  which  act 
to  alter,  or  perhaps  even  eliminate,  the  juncture  vortex  may  provide  fruitful.  Techniques 
which  act  to  reduce  or  eliminate  much  of  the  loading  asymmetry  near  the  root  by  altering 
the  effective  fin  attachment  angles  should  be  explored. 

7.2.2  Experimental.  It  is  recommended  that  the  experiment  at  Mach  4.9  be 
extended  to  include  a  hot-film  exploration  of  the  flowfleld  near  the  WAF  model.  The 
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data  from  such  experiments  would  provide  the  same  type  of  detailed  flow  angularity  and 
turbulence  data  that  were  obtained  at  Mach  2.8.  The  AFIT  Mach  5  wind  tunnel  is  currently 
being  modified  to  enable  such  measurements. 

It  is  also  suggested  that  surface  flow  visualization  (oil  film  and/or  pressure- sensitive 
paint)  be  obtained  at  several  Mach  numbers  on  both  spinning  and  non-spinning  WAF 
configurations  to  identify  significant  changes  or  bifurcations  in  the  flow  structure  near 
the  juncture.  Such  flow  visualizations  may  also  provide  insight  as  to  the  effects  of  unit 
Reynolds  number  on  the  flow  structure. 

7.2.3  Numerical.  From  a  numerical  prospective,  a  reasonable  next-step  would 
be  to  study  the  effects  of  cross-flow  on  the  flowfield  near  the  fin.  WAFs  most  certainly 
encounter  cross-flow  due  to  missile  spinning  and  yaw  (among  other  factors),  which  are  likely 
to  have  significant  effects  on  fin  loading,  perhaps  further  increasing  the  predicted  effect  of 
Mach  number  on  the  rolling  moment.  Experiments  on  sharp-fin/flat-plate  geometries  have 
demonstrated  a  considerable  dependency  of  juncture  vortex  structure  on  cross-flow.  I33'  34] 
The  effect  of  cross-flow  could  be  numerically  investigated  for  the  current  geometry  by 
solving  the  governing  equations  subject  to  periodic  boundary  conditions. 

Clearly  numerical  simulations  which  address  fin  interaction  issues  as  well  as  the  effects 
of  Mach  number  and  Reynolds  number  should  be  conducted.  Further,  simulations  which 
employ  turbulence  models  capable  of  resolving  flowfields  which  are  characterized  by  regions 
of  large  pressure  gradient  should  be  evaluated  against  the  current  body  of  experimental 
data. 
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Appendix  A.  Turbulent  Navier-Stokes  Equations 

A.l  Overview 

It  has  become  clear  that  computational  fluid  dynamics  (CFD)  will  play  an  increasing 
role  in  the  design  of  future  high  speed  vehicles  and  weapons.  In  fact,  due  to  the  difficulties 
and  large  expense  associated  with  ground  testing  such  configurations,  CFD  is  becoming  a 
necessity  in  the  design  process.  The  flow  over  such  vehicles  is  characterized  by  turbulent 
structures  that  are  not  well  understood.  This  lack  of  understanding  has  remained  the 
major  obstruction  to  accurately  simulating  complicated  high-speed  flows,  and  turbulence 
modeling  has  remained  the  controlling  factor  in  the  accuracy  of  predicting  such  flows.  Re¬ 
cently,  Settles  and  Dodson  l96l  pointed  out  that  the  modeling  community  is  sorely  lacking 
adequate  experimental  data  that  meets  the  high  standards  required  by  modern  code  val¬ 
idation.  One  of  the  objectives  of  the  experimental  portion  of  the  research  was  to  help  fill 
that  void. 

Direct  simulations  (of  relatively  low  Reynolds  number  flows  over  simple  geometries) 
have  shown  that  the  Navier-Stokes  equations  are  indeed  the  governing  equations  for  turbu¬ 
lent  flows.  However,  since  turbulent  flows  are  characterized  by  temporal  and  spatial  scales 
that  range  over  several  orders  of  magnitude,  direct  numerical  solutions  of  the  unsteady 
Navier-Stokes  equations  for  turbulent  problems  of  practical  interest  are  highly  unlikely  in 
the  near  future.  Thus,  researchers  are  restricted  to  considering  time  averages  of  turbu¬ 
lent  motion,  and  to  rely  on  approximation  methods  to  provide  solutions  to  high  Reynolds 
number  problems. 

The  two  most  prevalent  forms  of  the  governing  equations  are  the  Reynolds  (time) 
averaged  and  Favre  (mass- weighted- time)  averaged  forms  of  the  N-S  equations  (RANS  and 
FANS,  respectively).  In  either  case,  additional  fluctuation  cross- correlation  terms  appear 
in  the  averaged  form  of  the  equations,  resulting  in  a  mathematical  system  having  more 
unknowns  than  equations  (details  follow  this  section).  Reducing  the  number  of  unknowns 
to  the  number  of  equations  is  known  as  the  “closure”  problem. [68]  Thus,  the  function  of 
turbulence  modeling  is  to  accurately  represent  for  these  terms,  either  by  expressing  them 
as  functions  of  mean  flow  properties,  or  by  expressing  them  in  terms  of  additional  transport 
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equations.  The  form  of  the  closure  model  depends  on  the  type  of  averaging,  Favre  (FANS) 
or  Reynolds  (RANS).  In  the  past,  researchers  have  often  made  use  of  ad  hoc  assumptions 
to  achieve  closure  of  the  N~S  equations.  The  assumptions  result  in  turbulence  models 
that  are  rigorously  incorrect.  While  these  simplified  turbulence  models  do  serve  a  purpose 
and  have  had  some  degree  of  success,  they  are  inadequate  for  predictions  of  compressible 
viscous/inviscid  interactions  such  as  shock/boundary-layer  interaction.  t66>  961 

The  Reynolds  averaged  equations  are  obtained  by  un-coupling  the  instantaneous  flow 
properties  into  a  time- average  mean  value  plus  a  fluctuating  turbulent  contribution,  e.g., 

0  =  0+4/ 


where  0  is  the  time-averaged  quantity,  and  0'  is  the  instantaneous  turbulent  fluctuat¬ 
ing  component.  The  Favre- averaged  quantity  and  the  Favre  turbulent  fluctuation  of  the 
quantity  are  given  by 

0  =  ^  and  0"  =  0  —  0 
P 

respectively.  Note  that,  by  definition,  the  time  average  Reynolds  fluctuation  0',  is  zero, 
however,  the  time  average  of  the  Favre  fluctuation,  0",  is  non-zero.  The  form  of  the 
compressible  FANS  cross- correlation  terms  are  very  similar  in  appearance  to  those  of  the 
incompressible  RANS.  This  coupled  with  Morkovin’s  hypothesis^  which  states  that  “the 
turbulence  structure  is  unaffected  by  compressibility  as  long  as  the  fluctuation  Mach  num¬ 
ber  is  less  that  unity”  has  led  to  the  current  virtually  universal  trend  of  adopting  the 
FANS  equations  for  high  speed  compressible  flows.  As  a  result,  practically  all  compress¬ 
ible  turbulence  models  represent  direct  extensions  of  incompressible  formulations,  where 
the  constants  are  adjusted  and  the  density  is  allowed  to  vary.  Such  corrections  do  little 
more  than  correlate  the  data  upon  which  the  models  are  based. I75!  In  fact,  Morkovin’s  own 
Mach  1.77  expansion  fan/boundary  layer  interaction  data^  suggests  that  the  compressible 
upfvf  term  in  the  Reynolds  shear  stress  is  of  the  same  order  as  the  typical  incompressible 
pu'vf  term. 
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A. 2  Navier- Stokes  Equations 

In  Cartesian  coordinates  the  Navier-Stokes  equations  may  be  written  t6^: 
Conservation  of  Mass 

P,t  +  (puj),j  =  0 

Conservation  of  Momentum 

(pUi),t  +  (pUiUj),!  =  oijtj 

Conservation  of  Energy 

(pE),t  +  (pEui)ti  =  ((TijUj)'i  -  qiti 

or 

(pH),t  +  (pHui),i  =  P,t  +  (TijUj),i  -  qiti 

where 

E  =  e  +  \uiU{  ,  H  =  e  +  %  +  §«,•«,•  ,  e  =  CVT 

(Tij  —  P^ij  "l"  Tij  ,  T{j  —  2/J.Sij  fl  UlcfiSij  ,  S{j  —  2  (^i,j  "t"  ^J,i) 
p*  ~  \p  pbuik  j  q%  — 

A.S  Reynolds  (Time)  Averaged  Navier-Stokes  Equations 

When  using  Reynolds  averaging,  any  given  primitive  flow  quantity,  <j>,  is  expressed  as 
the  sum  of  a  time-averaged  component,  <f>  and  a  fluctuating  component,  <j>',  as  follows!6^: 

p  =  p  +  p'  Ui  =  Ui  +  u'i  p  =  p  +  p'  T  =  T  +  T1  E  =  E  +  E1 
H  =  H  +  H’  p  =  p  +  p'  p*=jF+p*'  CV  =  C^+C'V  Cp  =  ^+C; 

The  time  average  is  defined  by 

_  1  fto+T 

<£=-/  4>dt  (A.6) 

1  Jto 

where  the  characteristic  time  T  is  long  in  comparison  to  the  cycle  times  of  the  fluctuating 
component,  but  still  permits  gradual  time-dependent  fluid  motion  for  non-stationary  flows. 


(A.1) 

(A.2) 

(A.3) 

(A.4) 

(A.5) 


A- 3 


To  obtain  the  mean  flow  equations,  these  expressions  for  the  decomposed  variables 
are  substituted  into  the  governing  equations.  These  instantaneous  equations  are  then  time 
averaged,  expanded,  and  simplified. I105- 1161  The  results  are  the  time  averaged  (or  Reynolds 
averaged)  Navier-Stokes  equations  given  below  in  conservative  form. I38’  66'  131 
Mean  Continuity 

P,t  +  (pUj~™J),j  =  0  (A. 7) 

Mean  Momentum 

(P  W  +  P'u'i),t  +  (pW.uJ)j  =  ~P,i  +  ( Tij  +  rj  )j  (A.8) 

Mean  Energy 

(peo  +  P%),t  +  {phouj),j  =  (uirij  +  u'Tij  ~  9i  “  ?J),j  (A-9) 

with  the  Equation  of  State 

p  =  pRT  +  Rp T7  (A. 10) 

where  the  compressible  RAN S  turbulence  terms  have  been  defined  to  be 

mj  =  -p'u'j 

Tjj  =  —p  u'iu'j  -  ul  p'u'j  —  uj  p'u'i  -  p'u'iu'j  (A-11) 

qj  =  +p  h'0u'i  +  h0  p'u'i  +  u~i  P'K  +  p'h'Qu'i 

However,  this  process  introduces  an  additional  6  unknowns.  In  these  equations  e0  and  h0 
are  stagnation  conditions  (i.e.,  e0  =  e  +  |u,-u,).  The  terms  which  are  due  to  Reynolds 
averaging  are  the  turbulent  apparent  mass  flux,  mf ,  compressible  turbulent  shear  stress, 
r,”,  and  compressible  turbulent  heat  flux,  qj .  For  thin  layer  flows,  these  reduce  to: 

rriy  =  -pV  rjv  =  -p  vFv'  -urnFy  qy  =  +p  h'Qv'  -  rn?y  (A. 12) 
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and  for  incompressible  flows: 


mj  =  0  =  -p  u'iu'j  qj  =  +p  h'0u'i  (A.13) 


A. 4  Favre  Averaged  Navier-Stokes  Equations 

When  employing  Favre  averaging,  any  given  primitive  flow  quantity,  4>,  is  expressed 
as  the  sum  of  a  mass-weighted-averaged  mean  value  (Favre  averaged),  <f>,  and  a  fluctuating 
component,  <f>".  The  Favre  averaged  quantity  and  the  Favre  turbulent  fluctuation  of  the 
quantity  are  defined  as  I66l : 

4>  =  ^  and  <j>"  =  <f>  -  (A.14) 

P 

It  should  be  noted  that  although  <j>"  ^  0,  it  is  also  true  that 

~p<i>  =  p{4>  +  4>")  =  p4>  +  p4>"  =  ~p4>  +  p4>"  =  p^=r  +  p<t>"  -  p<t>  +  p4>" 

so,  it  must  be  true  that  p<j>"  =  0.  Also,  the  following  relationships  can  be  derived 


-  -  „  ,  —  p'4>" 

$  -  d>  =  d>"  -  4>'  =  <t>"  =  -~- 

p 


p 


The  Favre  averaged  equations,  as  presented  by  Marvin  I68?  (and  corrected  here),  axe 


Continuity: 


dp  d(puj)  _ 
8t  +  dxi 


(A.15) 


Momentum:  _ 

d{puj)  djpujUj)  _  dp  d(Ttj  -  pu-u'-) 

dt  dxj  dxi  dxj 

Energy: 


(A.16) 


dph  d[pMj)  _  dp  .  djp_  „  dp  d(-q}  -  pti'u")  Qu" 
dt  dxj  dt  dxj  Ui  dxj  dxj  T,J  dxj 


(A.17) 
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These  equations  have  the  same  form  as  the  incompressible  Reynolds  averaged  equa¬ 
tions,  except  that  the  Reynolds  stresses  pu"u"-  include  the  density  fluctuations,  which 
must  be  accounted  for  in  some  way.  Compressibility  is  generally  included  by  replacing  the 
density  with  the  mean  density  neglecting  the  terms  which  involve  additional  correlations 
arising  from  compressibility.  Researchers  using  the  Favre  forms  of  the  governing  equations 
cite  that  avoiding  some  of  the  fluctuating  density-generated  terms  simplifies  closure  of  the 
system.  However  since  much  of  the  experimental  data  provides  RANS-type  information, 
modeling  the  effects  of  compressibility  is  problematic.  The  lack  of  an  explicit  fluctuating 
density  component  in  the  turbulent  shear  stresses  is  one  of  the  specific  shortcomings  of 
this  model. 

The  compressible  FANS  turbulence  terms  may  be  expressed  as 

mI  =  0  T J  =  -pu'fu'!  qj  =  ph'lu'l  (A. 18) 

Note  that  these  compressible  Favre  averaged  terms  are  similar  in  appearance  to  their 
incompressible  Reynolds  averaged  counterparts  (Equation  A. 13).  Bowersox  and  Schetzt16) 
have  shown  that  cross- wire  anemometry  is  well  suited  for  measurement  of  RANS  turbulent 
terms.  In  particular,  they  have  shown  that  the  total  shear  stress  can  be  directly  measured 
for  thin  layer  type  flows.  Furthermore,  if  the  effects  of  the  pressure  fluctuations  on  the  hot¬ 
wire  response  are  small,  then  the  multiple  overheat  cross-wire  results  can  be  decomposed 
into  all  of  the  terms  in  Equation  A.ll.  This  assumption  has  been  verified  by  Bowersox 
and  Schetzt16)  for  a  Mach  4.0  free  mixing  layer,  and  Kistlert53)  suggests  that  is  valid  for 
supersonic  boundary  layers  up  to  Mach  4.7. 

There  are  also  forms  of  the  Navier- Stokes  equations  in  which  both  Reynolds-  and 
Favre  averaging  are  used!66)  In  these  equations,  the  Reynolds  average  is  often  used  for  the 
density,  pressure,  transport  coefficients,  and  specific  heats,  and  Favre  averaging  is  used  for 
the  other  quantities. 
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A. 5  Zero-Equation  Turbulence  Models  (Baldwin-Lomax) 

Both  the  RANS  and  FANS  equations  have  additional  fluctuation  cross-correlation 
terms  appear  in  the  averaged  form  of  the  equations,  causing  them  to  have  more  unknowns 
than  equations.  Thus,  turbulence  modeling  must  be  used  to  reduce  the  number  of  un¬ 
knowns  (or  increase  the  number  of  equations)  to  close  the  system.  This  can  be  accom¬ 
plished  by  expressing  the  cross- correlation  terms  as  functions  of  mean  flow  properties,  or 
by  expressing  them  in  terms  of  additional  transport  equations. 

The  following  discussion  is  in  no  way  intended  to  be  a  thorough  review  of  turbulence 
modeling,  but  only  a  brief  overview  of  the  type  of  turbulence  model  used  in  the  current 
research. 

The  simplest,  and  most  widely  used  turbulence  models  are  derivatives  of  the  method 
of  Cebeci  and  Smith  l18’ 19'  17L  The  main  drawback  of  this  method  is  that  the  boundary 
layer  thickness  has  to  be  known.  In  zero-equation  (algebraic)  models,  the  concept  of 
a  turbulent  or  eddy  viscosity,  fit,  is  used.  This  eddy  viscosity  is  simply  added  to  the 
molecular  viscosity  in  the  governing  equations,  i.e., 

fi  -►  fi  +  fit 

Also,  a  turbulent  Prandtl  number,  Prt,  is  defined  such  that 

JL  t  El 

Cp  Pr  +  Prt 

in  the  energy  equation.  Baldwin  and  Lomax!11)  modified  the  model  to  eliminate  the  need 
for  finding  the  edge  of  the  boundary  layer,  S. 

f 

(Mt)inner  V  ^  Vc 

f*t  = 

(^)outer  Vc  <  V 

where  y  is  the  normal  distance  from  the  wall,  and  yc  is  the  smallest  value  of  y  at  which 
the  outer  formulation  is  to  be  used. 
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The  inner  region  is  based  on  the  Prandtl-van  Driest  formulation 


(/pinner  — 

where  the  mixing  length  and  scaled  coordinate  are 

t=ky\  1  -  e~y+'A+}  and  y+  = 

L  J  /it v  Hw 

where  ur=^/^  is  the  friction  velocity,  and  |w|  is  the  magnitude  of  the  vorticity.  The  outer 
layer  model  is  a  modification  of  the  Clauser  formulation. 

(/^f  ) outer  —  K  Cep  p  Fyr& Ice  Tkleb(2/) 

in  which  K  is  the  Clauser  constant,  Cep  is  a  constant,  and  the  wake  function  is  given  by 


and  j/max  is  the  value  of  y  at  which  jFmax  occurs.  The  KlebanofF  intermittency  factor  is 
given  by 

Tideb(y)  =  1  +  5.5 

\  Vma.x  ) 

The  constants  specified  to  agree  with  Cebeci-Smith  formulation  for  constant  pressure 
boundary  layers  at  transonic  speeds,  and  are  given  as 


A+  =  26  Cep  =  1.6  Ckieb  =  0.3  Cwk  =  0.25 
k  =  0.4  K  =  0.0168  Pr  =  0.72  Prt  =  0.9 
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While  simulating  the  flow  about  supersonic  pointed  bodies  at  large  incidence,  Degani 
and  Schiff  modified  the  method  to  properly  account  for  the  associated  large  regions  of  cross- 
flow  separation.  t22l  This  allowed  the  method  to  more  accurately  predict  these  turbulent 
3-D  vortical  flows  using  the  thin  layer  parabolized  Navier-Stokes  (PNS)  equations. 
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Appendix  B.  Hot-Film  Methods  for  Turbulence 

Care  must  be  taken  when  comparisons  between  numerical  and  experimental  results 
are  made.  However,  it  is  clear  that  hot-film  anemometry  responds  to  time-averaged 
mass  fluxes  and  total  temperatures!14].  The  present  research  uses  single-overheat  ther¬ 
mal  anemometry  to  measure  the  compressible  Reynolds  turbulence  data  in  a  Mach  2.8 
(Re/£=18x  106m-1)  flow  involving  a  shock/boundary-layer  interaction.  Detailed  three- 
dimensional  surveys  of  pressure  and  mass  fluxes  have  been  obtained,  enabling  the  calcu¬ 
lation  of  turbulent  quantities  such  as  the  turbulence  intensities  and  Reynolds  turbulent 
shear  stresses. 

The  constant  temperature  hot-film  anemometer  I84'  671  records  the  voltage  required  to 
maintain  the  film  at  a  constant  known  temperature.  The  power  required  to  maintain  this 
temperature  is  equivalent  to  the  heat  transfer,  qj,  between  the  hot-film  and  the  surrounding 
flow.  The  Nusselt  number  can  be  related  to  the  heat  transfer  from  the  film  by 


Nu  = 


U 

irk£f(Tf  —  Te) 


(B-l) 


Where  Te  is  the  temperature  the  unheated  film  would  approach  under  these  specific  flow 
conditions  (equilibrium  temperature). 


The  basic  shape  of  a  hot  film  probe  used  in  this  research  is  a  cylinder  and  the  form 
of  the  Nusselt  number  for  compressible  flow  in  dimensionless  heat  transfer  is  t57]; 


Nu  =  fcn(£f/d,  M,Pr,Ree,r) 


(B.2) 


lf/d  is  the  film  aspect  ratio;  M  is  the  Mach  number;  Pr  is  the  Prandtl  Number;  Ree  is  the 
effective  cooling  Reynolds  number  base  on  film  diameter;  and  r  is  the  temperature  loading 
factor.  Th e  temperature  loading  factor  can  be  expressed  as 


T  = 


Tf-Tt 

Tt 


where  Tf  is  the  film  temperature  and  Te  is  the  equilibrium  temperature,  or  the  temperature 
that  the  unheated  film  would  attain  if  placed  in  the  flow.  For  Reynolds  numbers  greater 
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Figure  B.l  Sample  two- film  calibration  using  King’s  Law. 

than  about  20,  Te  is  about  97%  of  T*.  When  the  Mach  number  normal  to  the  film  is 
greater  than  1.2,  or  M simp  >1,  Pr  is  assumed  constant,  and  the  aspect  ratio  lj/d  1, 
the  function  for  the  Nusselt  number  simplifies  to 

Nu  = /cn(Ree,r)  (B.3) 

It  has  been  determined  experimentally  by  Bowersox  and  Schetzl16)  that  King’s  Law,  the 
functional  relationship  between  Nusselt  number  and  Reynolds  number  for  incompressible 
flow,  is  also  an  acceptable  relationship  for  compressible  flows. 

Nu  =  akVfc^+bk  (B.4) 


where  the  constants  a*  and  bk  must  be  experimentally  determined  for  each  value  of  T/. 
Sample  data  from  a  two-film  cross- wire  calibration  along  with  the  corresponding  King’s 
Law  curve  fits  are  presented  in  Figure  B.l.  By  definition,  the  Nusselt  number  is  also 
proportional  to  the  film  power 


Nu  = 


U 

7 rktef(Tf  -  Tt) 


(B.5) 
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where 


jp 

q,  =  i?R,  and  i,  =  +  + 

Here  Rf  is  the  film  resistance,  Rs  is  the  resistance  of  the  resistor  in  series  with  the  film, 
and  Rl  is  the  probe  lead  resistance.  Assuming  that  Te=Tt  (which  results  in  minimal  error 
if  done  in  both  the  calibration  and  the  data  reduction),  the  Nusselt  number  can  be  written 
as 

■KT  _  _ •£*/  _ 1 _  /T>  g\ 

1NU  "  (Rf  +  RS  +  Rtf  i rkt£f(Tf  -Te)  K  ’  ’ 

The  turbulent  power  laws  for  viscosity  and  thermal  conductivity  are  used  to  determine 
Nusselt  and  Reynolds  numbers. 


kt  =  nk  =  0.89 

Pt  =  (§)""  »/•  =  O-7^ 

Combining  (B.4),  (B.6)  ,  and  (B.7),  gives  the  hot  film,  response  equation 

'Tt\-n»/ 2 


E/  =  (T1\ 
Co  \T0J 


nk 


afc^Re0e 


+  bk 


(Tf  —  Tt) 


where  Re0e  is  the  effective  Reynolds  number  with  p  =  po,  and 

(Rf  +  Rs  +  RLf 


C0  = 


Rf 


-irlfk o 


(B.7) 


(B.8) 
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B.l  Turbulence  Fluctuations 

Replacing  Ef,  Reoe,  and  Tt  by  their  mean  and  fluctuating  components  in  the  hot 
film  response  equation  (B.10),  using  the  Binomial  Theorem,  retaining  only  the  first  order 
terms,  and  noting  that 


eI  = 

Co 


*  [ajfe-v/5^7 


+  bk 


(Tf  —  Tt) 


(B.10) 
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then  solving  for  Ej/Ef7  the  hot  film  fluctuation  equation  is  derived  in  the  form 


where  the  hot  film  sensitivities  are  given  by 


1+ 


bk 


Q'kyR'^e  J 


-1 


— Tt  njs 

and  g  =  — - =-  H - /  na 

3  c\f  m  m  \  ■  <\  j  V- 


2  (Tf-T,)  2 


(B.ll) 


(B.12) 


For  the  single-overheat  technique  (141 ,  the  sensitivity  to  total  temperature  fluctuation  is 
minimized  by  operating  the  hot-film  at  large  overheat  ratios  «2);  hence  g  is  small. 
For  flows  where  the  total  temperature  fluctuations,  T/,  is  also  small,  the  second  term 
in  Equation  B.ll  can  be  neglected.  For  Mach  3  boundary  layers,  the  total  temperature 
fluctuation  has  been  found  to  be  small  (^«2-3%),  thus  neglecting  that  term  is  reasonable. 
Specifically,  the  AFIT  Mach  3  wind  tunnel  has  been  found  to  maintain  total  temperature 
fluctuations  below  about  2.0%J44'  45'  70’  7S!  To  evaluate  /,  Ree  is  found  from  Equation  B.10, 
where  T%  is  assumed  equal  to  the  plenum  total  temperature,  Ttx . 


B. 1.0.1  Decomposition  into  Cartesian  Coordinates.  To  derive  the  formulas 
for  analysis,  the  effective  Reynolds  number  given  by  the  previous  section  must  be  related 
to  the  x-y  coordinate  system.  Since  <p  did  not  vary  substantially  from  the  calibration 
values,  the  cosine  law  was  assumed  valid,  i.e., 


Ren 

^  j 

►  — 

COS<£> 

sing? 

Rex  1 

Ret  J 

—  sin  9? 

cos  cp 

Rey  J 

Thus  the  effective  cooling  Reynolds  number  becomes 


(B.13) 


R-6e2  —  AiRex2  +  2A2RexRey  -f  A^Rey2 


(B.14) 


B-4 


where  A*  are  given  by 


A\  —  cos2  cp 

A2  =  cosy?  sin  <p  (B.15) 

A3  =  sin2  <p 

and  (f  is  the  incidence  angle  of  the  hot  film  to  the  flow. 

Replacing  Reoe ,  Reo,  and  Re0j,  by  their  mean  and  fluctuating  components  and  ap¬ 
plying  the  binomial  theorem,  with  Rq  <  1,  one  can  show  that 

Reo,.,  =  Re0:t  \J Bzj  (B.16) 


and  write  the  fluctuation  equation  as 


.  „  /Be o,'\ 


(B.17) 


where  the  index  j  sums  over  the  two  films  on  the  cross-film  probe,  and 


Rq  = 
Bi  = 
B2  = 
Bz  = 


Al 

b3 

62. 

B3 

A\  +  2A2Rq 


(B.18) 


Solving  this  set  of  equations  and  decomposing  into  x  and  y  components 


Re0l  = 


•2  Reoel  /A21  —  Reoe2  /A22 


A11/A21  —  A12/A22 


(B.19) 


Reo..  = 


_  1  Repel  Mu  —  Rep^,  ! A\2 

2Re0l  A21/A11  -  A22/A12 
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The  turbulence  variables  can  be  decomposed  into  x  and  y  components  via  (B.17)  as 


WV  1 

Re0l  )  D22 

Reo/\2_  1 
Re0l )  Di2 


1  ( Re0/\2  2  /Rep/\  StepA 

B2\ 2  \  Reoe  /  a  B2\B22  y  Re0(.  ) l  \  Re0e  /  2  B222  \  Re0e  /  2 

_  _ (B.20) 

1  /Reo/\2  2  7y\  fRV^  1  /Reo/V 

Bn2V»®o./i  ^n^VReoJi  UeoJ2  5122^Re0J2 


1 

2B\\B2i 


(B.21) 

(B.22) 


where  D\  =  ( B2i/Bn  -  B22/Bi2)  and  D2  =  ( Bu/B2i  -  Bi2/B22).  The  covariance  of  the 
two  films  can  be  expressed  as 


(B.23) 


For  the  x-z  plane,  the  above  equations  are  used  with  w  replacing  v,  and  z  replacing  y.  In 
this  case,  the  cross-film  is  rotated  90  degrees  (i.e.  an  x-z  cross-film  is  used). 


B.2  Turbulence  Transformation  -  Reynolds  Shear  Stress 

Bowersox  and  Schetzl16!  have  shown  that  the  Reynolds  shear  stress  can  be  expressed 
in  terms  of  the  conservative  cross-film  variables  as 


(puiYjpujY  | 

P 


“’">(0 


where  the  second  term  on  the  right  hand  side  has  been  shown  to  be  much  less  than  the  first 
term  for  thin  layer  type  flows.  For  the  present  flowfield,  the  flow  angle  was  usually  less  than 
5-10°,  hence  neglecting  the  second  term  would  have  probably  been  acceptable.  In  fact, 
the  second  term  was  calculated  using  the  methods  of  Bowersoxt1^  and  was  determined  to 
be  at  least  an  order  of  magnitude  smaller  than  the  first  term  for  all  of  the  data  presented. 
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Appendix  C.  Experimental  Uncertainty  Analysis 


Error  and  uncertainty  are  inherent  to  experimental  research.  This  section  attempts  to 
identify  the  possible  sources  of  measurement  error  in  the  experiments,  and  then  quantify 
the  effects  of  these  errors  on  the  presented  data.  Based  on  the  analysis  presented  by 
Bowersoxt1^  the  Euclidean  (Z2)  norm  is  utilized  to  assess  the  cumulative  effects  of  error 
sources.  The  L\  norm  (a  summation  of  the  absolute  error  values)  has  been  found  to  be  too 
conservative  an  estimate  when  compared  to  a  perturbation  analysis  of  the  data  reduction 
equations.  The  i2  norm  of  a  set  data,  ®,-,  is  defined  byt102! 


IMIz  =  11*1, *2, 


Xnh  = 


Thus,  the  total  dimensional  error  is  defined  to  bet42! 


(C.l) 


(C.2) 


where  the  index  i  runs  over  the  various  measurement  errors,  ZXi,  associated  with  determi¬ 
nation  of  the  result,  R.  The  dimensional  errors  can  be  normalized  by  reference  values  to 
obtain  nondimensional  errors  (i.e.  Sr  =  £r/R).  The  effects  of  measurement  errors  were 
propagated  through  the  data  reduction  process,  where  the  equations  were  linearized  to 
provide  approximate  error  bounds  on  the  processed  data  for  the  experiments. 


C.l  Measurement  Errors 

Every  measurement  has  an  associated  error.  It  is  assumed  all  measurement  errors 
are  random,  with  a  Gaussian  distribution.  The  assumption  of  random  errors  precludes  the 
existence  of  biased  errors  or  blunders. 

C.1.1  Conventional  Probes.  Pressure  was  measured  by  transducers  which  mea¬ 
sure  gauge  pressure.  An  error  of  ±0.0034atm  was  assumed  in  measurement  of  ambient 
pressure,  determined  by  the  smallest  increment  on  the  barometer. 
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The  pressure  transducers  used  to  measure  pitot,  cone-static  and  plenum  pressure  are 
advertised  to  have  accuracies  of  0.5%,  0.5%  and  0.4%,  respectively.  t3°l  Calibration  and 
digital  conversion  of  the  pressure  data  for  storage  and  processing  also  adds  error.  All 
pressure  readings  were  conditioned  the  same  style  of  indicator  (Endevco  4428A)  with  a 
maximum  gain  error  of  ±0.5%  and  gain  stability  of  ±0.2%.l31l  The  probe  pressure  units 
were  re-zeroed  daily,  and  never  deviated  more  than  O.Olpsi  or  O.lpsi  for  the  probe  and 
plenum  pressures,  respectively.  Finally,  Volluzt115'  reports  that  turbulence  induces  about 
±0.0068atm  error  for  both  pitot  and  cone  static  probe  types. 

Using  multiple  overheat  (MOH)  hot-film  anemometry,  the  AFIT  Mach  3  wind  tun¬ 
nel  has  been  found  to  have  boundary  layer  total  temperature  root  mean  square  (RMS) 
fluctuations  of  2.0%J74'  691  The  RMS  fluctuations  are  included  as  errors  in  this  analysis  as 
they  provide  a  definitive  bound  of  measurement  uncertainty  from  large  data  samples.  The 
plenum  total  temperature  was  observed  to  vary  no  more  than  ±2.36K  during  tunnel  runs. 
The  plenum  thermocouple  and  display  unit  were  accurate  to  within  ±1K. 

The  degree  of  flexing  experienced  by  each  type  of  probe  during  tunnel  operation  was 
determined  by  aligning  two  grids  onto  the  Plexiglas  plates  to  measure  x  and  y  position. 
The  flex  angle  and  position  were  then  determined  as  a  function  of  initial  measured  position, 
and  accurate  to  within  1°  (1.1mm).  The  measurement  of  x  and  y  position  was  accurate 
to  ±0.5mm.  Without  Plexiglas  windows  on  the  top  and  bottom  of  the  tunnel,  it  was 
more  difficult  to  determine  the  z  position  of  each  probe.  The  error  in  z  position  was 
assumed  to  be  twice  the  measurement  error  in  x  or  y  (1.0mm).  An  additional  error  from 
digitization  occurred  in  the  measurement  of  y  position  -  as  y  position  was  recorded  using  a 
linear  displacement  voltage  transducer,  LDVT.  Since  the  LDVT  was  sampled  with  a  12.0V 
range,  the  above  error  in  digitization  voltage  was  multiplied  by  the  calibration  slope.  This 
manipulation  resulted  in  an  error  in  y  position  from  digitization  of  ±0. 015mm.  Errors  in 
position  were  also  incurred  due  to  the  flexing  of  the  probes  during  tunnel  operation. 

C.1.2  Hot-Film  Probes.  The  possibility  of  measurement  error  is  greatest  in  the 
angular  rotation  of  the  probe.  Since  the  probes  are  calibrated  at  an  upstream  location 
then  moved  to  the  downstream  measurement  locations,  great  efforts  were  taken  to  ensure 
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that  the  probe  remained  aligned  with  the  tunnel  axis.  Probe  alignment  was  “spot  checked” 
where  possible  by  examining  the  measured  flow  angularity  in  the  freestream. 

Hot-wire  measurements  have  shown  1.0%  root  mean  square  fluctuation  in  voltage 
for  freestream  flow.  Additionally,  a  digitization  error  of  ±0.003V  was  accounted  for  in 
the  12.0V  sampling  range.  The  errors  in  pressure,  temperature,  position  and  voltage  have 
been  summarized  in  Table  C.l. 


Table  C.l  Measurement  error  bounds. 


Measured  Property,  x 

Error 

Mach  2.8 

Mach  4.9 

X 

1.2mm 

1.2mm 

y 

0.5mm 

0.5mm 

Z 

1.0mm 

1.0mm 

Tn 

2.3% 

2.2% 

Pti 

0.8% 

0.7% 

Pt2 

0.9% 

0.8% 

P zs 

3.5% 

1.3% 

Ef 

1.0% 

— 

C.2  Error  Propagation 

The  errors  listed  in  Table  C.l  have  an  influence  on  all  subsequent  data  reduction.  To 
determine  the  influence  on  calculations,  the  equations  used  are  linearized  about  freestream 
conditions.!44!  Tables  C.2  and  C.3  list  the  freestream  conditions  and  typical  freestream  hot- 
film  calibration  parameters  used  in  the  analysis,  as  well  as  to  normalize  the  uncertainty 
estimates. 

C.2.1  Properties  Determined  with  Pressure  Probes.  The  majority  of  mean  flow 
calculations  are  based  on  the  local  Mach  number,  where  the  Mach  number  is  calculated 
from  a  curve  fit  to  experimental  data  (Equations  3.1  and  3.2).  As  an  example  of  error 
propagation,  the  analysis  is  applied  to  Equation  3.1  which  is  of  the  form 

I  =  Cb  +  +  C2e  +  C3|3  +  C4|4  =>  M  =  M2[C0  +  Cal  +  C2|2  +  C3|3  +  C4|4] 
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Table  C.2  Freestream  conditions. 


Condition 

M=2.8  Value 

M=4.9  Value 

Tn 

294K 

354K 

Pti 

2.085atm 

32.0atm 

to 

0.7175atm 

1.86atm 

Pcs 

bhkkqq 

0.3388atm 

a 

mmm 

0.1725 

M 

2.80 

4.90 

7 

1.4 

1.4 

Rair 

■KME3 

T 

114.5K 

61.5K 

aoo 

■m 

MBSBEM 

Pi 

0.065atm 

0.070atm 

pi 

U 

■H 

766? 

P\U 

Reoo 

1.8xl06/m 

Table  C.3  Typical  freestream  hot-film  parameters  (M=2.8). 


Condition 

Value 

Ef 

5.0V 

Tf 

700.0K 

Te 

294.0K 

Rf 

10.0ft 

Rs 

50.0ft 

Rl 

0.0ft 

a,k 

0.11 

bk 

d 

5.1  •  10~5m 

L 

0.001m 

Ree 

300 

Nu(wire) 

1.8 
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Application  of  Equation  C.2  results  in  the  normalized  error  at  Mach  2.8  of 


£m  =  ^  =  ||M(Cif  +  2 C2e  +  3Cse3  +  4 C4t4)  •  e€||2  =  0.725s*  (C.3) 

Similarly,  Equation  3.2  at  Mach  4.9  gives 


eM  =  ||M(Cif  +  C2I)  -£$||2  =  0.855s*  (C.4) 


The  pressure  ratio,  £  is  the  ratio  of  pressures  measured  with  a  10°  or  20°  cone- static 
probe  and  a  pitot  probe. 

p  i  p.  i 

(C.5) 


P _  P :5 

*~pTi 


. Ptl 

cs  Pt2 


t2 


Two  different  plenum  pressures  (Pn\t2  and  Pn\cs)  are  used  to  minimize  errors  caused  by 
differences  between  the  two  separate  tunnel  runs  (one  for  each  probe).  Since  f  was  a 
combination  of  pressure  measurements,  the  errors  combine  with  the  Euclidean  norm 


H  =  p£Pni£Pt2,£pJ\2  (C-6) 

The  Mach  number  and  plenum  temperature  were  used  to  derive  the  local  tunnel 
temperature,  T\,  assuming  isentropic  flow.  Although  isentropic  flow  was  violated  by  shocks, 
no  method  for  measuring  local  temperature  was  available.  The  equation  of  state  and 
ideal  gas  laws  were  then  used  to  calculate  density  and  speed.  Table  C.4  summarizes  the 
propagation  errors  of  the  pressure  probe  measurements. 

C.2. 2  Properties  Determined  with  Hot- film  Measurements.  Single  overheat  (SOH) 
hot-film  anemometry  was  used  to  measure  mass-flux  mean  flow  and  RMS  fluctuations. 
SOH  analysis  assumes  negligible  total  temperature  fluctuations  and  has  proved  to  be  valid 
for  this  experimental  facility,  t74- 261  Additionally,  cross-wire  measurements  provided  flow 
direction  information.  A  detailed  explanation  of  the  hot-film  data  reduction  techniques 
which  were  used  is  presented  in  Appendix  B.  Due  to  the  complicated  nature  of  the  hot- 

film  data  reduction,  a  logarithmic/derivative  technique!44!  was  applied  to  the  hot-film  data 

* 

reduction  equations  to  estimate  the  propagation  of  errors  throughout  the  hot-film  analysis. 
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Table  C.4  Pressure  probe  related  error  bounds. 


Calculated 

Property 

Derivation 

Error 

H 

to 

bo 

o> 

II 

s 

M=2.8 

M=4.9 

£ 

1 1 2^-Ptl  5  &Pt2  ’  ^ Pc S  1 1 2 

BMijyMllIll 

3.9% 

2.0% 

M 

0.725^ 

2.8% 

1.7% 

Tx 

||£Xtl,1.22l£M||2 

||£Ttn  1.655£m||2 

4.1% 

3.6% 

a 

0.5£Xl 

O.SeTj 

2.1% 

1.8% 

Pi 

||£pt2,  1.907£m||2 

lkp,2,  1-970£m  b 

5.5% 

3.5% 

pi 

1  \ePnSTih 

lkfl,€li||2 

6.9% 

5.0% 

u 

||£m?£c*||2 

II^Mb  H2 

3.5% 

2.5% 

The  Nusselt  number  (of  the  hot-film),  Nu,  was  determined  from  the  power  consump¬ 
tion  required  to  maintain  a  constant  wire  temperature  (Equation  B.6).  The  error  in  Nu 
was  a  function  of  the  measured  voltage  error  and  total  temperature  error.  The  effective 
Reynolds  number,  Ree,  was  determined  by  a  curve  fit  of  Nu  ( King’s  Law ,  Equation  3.3). 
Additionally,  the  fluctuation  in  Ree  was  a  function  of  measured  voltage  error.  These  errors 
were  combined  under  the  L<i  norm.  The  error  in  mass-flux  is  equivalent  to  the  error  in 
Reynolds  number,  and  it  assumes  that  the  transverse  mass-fluxes  suffer  the  same  error  as 
the  axial  component.  Table  C.5  summarizes  the  propagation  errors  of  the  hot-film  probe 
measurements. 

C.2.3  Separation  of  Primitive  Fluctuations .  The  separation  equations  for  the 

primitive  fluctuations  were  linearized  about  the  reference  conditions  of  Table  C.2.  The 
results  of  variable  separation  error  analysis  are  summarized  in  Table  C.6. 


C.3  Comparison  of  u~v  and  u~w  Hot-film  Probes 

The  agreement  between  the  axial  mean  and  fluctuating  quantities  obtained  by  the 
two  probes  was  considered  excellent  (Figure  C.l).  The  small  differences  were  attributed 
primarily  to  the  high  flow  angles  (<£«±10°)  experienced  by  the  probes  which  have  finite 
(1mm)  wire  separations. 
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Table  C.5  Hot-film  related  error  bounds  (M=2.8). 


Calculated  Property 

Derivation 

Error 

Nu 

||2.0e£„erJ2' 

3.0% 

Ree 

6.1% 

/ 

0.5£Re<. 

3.1% 

9 

\^Ttl,£fh 

3.8% 

R ejf 

2.0£Ree 

12.2% 

Rex 

ll£Re2£1»£ReeiRee2,£Re2j|2 

21.2% 

Rg^ 

£Re2 

21.2% 

Re^Re* 

£Re2 

Rex 

°-^£Re2 

10.6% 

Rg^ 

10.6% 

Re* 

£Rex 

10.6% 

pu 

£Re, 

10.6% 

pv 

£Re„ 

10.6% 

pw 

£Re* 

10.6% 

Table  C.6  Turbulent  fluctuation  error  bounds  (M=2.8). 


Calculated  Property 

Derivation 

Error 

(pu)' 

£Re* 

10.6% 

(pu)' 

£Rev 

mmm 

(pw)' 

^Re* 

10.6% 

(pu)'(pv)' 

^Rcj  Rcy 

mm 

(pu)'(pw)' 

^RexRez 

21.2% 

p' 

||0.458£m,  0.77l£(^v||2 

8.3% 

u 1 

v' 

sut 

13.4% 

w' 

13.4% 
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[b]  Axial  Turbulence  Intensity,  ( pu)'/(pu ). 


Figure  C.l  Comparison  of  hot-wire  probe  orientation  at  x=0.69c  (M=2.8). 


C-4  Assessment  of  Probe  Location  Error  Using  Numerical  Solutions 

The  effects  of  probe  location  error  were  also  assessed  by  ‘probing’  the  computational 
results.  By  assuming  a  location  error,  then  examining  the  CFD  results  at  the  extremes 
of  that  error  band,  an  estimate  of  the  sensitivity  of  measured  flowfield  variables  to  probe 
location  error  may  be  ascertained.  Figure  C.2  shows  the  effect  that  a  ±lmm  probe  location 
error  in  the  x  and  2  directions  (in  the  numerical  solution)  has  on  the  measured  pitot 
pressure.  Due  to  the  high  transverse  pressure  gradients,  the  measured  data  are  more 
sensitive  to  variations  in  2  (with  measurement  errors  of  up  to  8.5%).  than  it  is  to  x  (up 
to  5.5%  on  the  concave  side).  In  reality,  it  is  likely  that  the  probe  experiences  less  flexing 
as  it  is  moved  away  from  the  model  wall.  If  this  is  true,  Ax  becomes  less  positive  as  Y/ 6 
increases,  and  a  profile  more  closely  resembling  the  experimentally  observed  data  (without 
the  acceleration  peak  over  Y/<5e[0.8,1.8])  could  result  on  the  concave  side  of  the  fin.  The 
measured  data  on  the  convex  side  would  be  relatively  unaffected  as  it  is  comparatively 
invariant  with  Ax. 

Many  of  the  differences  between  measured  and  predicted  flow  angularity  may  also  be 
attributed  to  probe  location  errors.  Figures  C.3  through  C.6  show  the  effects  of  a  ±lmm 
probe  location  error  in  the  x  and  2  directions  on  the  calculated  flow  angularities.  For  ex¬ 
ample,  at  a  location  only  1mm  closer  to  the  fin  than  estimated  on  the  convex  side,  the  nu- 
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Figure  C.2  Effect  of  a  variation  in  probe  position  on  computed  pitot  pressure  (M=2.8). 


merical  simulation  almost  exactly  reproduces  the  experimentally  observed  flattening  of  the 
horizontal  angularity,  6 ,  near  the  body  (0.3<Y/£o<0.8)  (c/Figures  5.8[c]  and  C.3[b]).  At 
this  same  location,  the  experimentally  observed  inflection  seen  in  the  azimuthal  flow  angle, 
<j>,  is  also  very  closely  reproduced  in  the  numerical  solution  (c/Figures  5.8[d]  and  C.4[b]). 

On  the  concave  side  of  the  fin,  moving  the  probe  location  1mm  further  away  from  the 
fin  (in  the  numerical  solution)  yields  the  drastic  reduction  in  the  horizontal  flow  angularity, 
0,  near  the  body  (Y/£o<0.7)  similar  to  that  observed  in  the  experiment  (c/Figures  5.9[c] 
and  C.5[b]).  The  sensitivity  of  6  to  probe  location  error  in  the  2  direction  in  this  region 
(z/czi—Al)  can  also  be  Figure  5.11[a].  Near  the  body,  the  flow  can  be  directed  either 
toward  or  away  from  the  body  with  a  small  variation  in  2. 
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Cv:  x=0.69c 
Cv:  x=0.69c+1mm 
Cv:  x=0.69c-1mm 


Cv:  Xs0.69c 

x=0.69c :  z=+0.47-1  mm 

x=0.69c :  zsr+0.47+1  mm 


[a]  Axial  Ax  =  ±lmm. 

[b]  Transverse  Az  =  ±lmm. 

Figure  C.3  Effect  of  a  variation  in  probe  position  on  0 ;  convex  side  (M=2.8). 

i  3.5 
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I  2*5 

i» 

; 1 

jjl 
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y  \+Vw  a—  v»o?v 

f  - Cv:  x=0.69c+1mm 

X  - Cv:x=0.69c-1mm 

A* - x=0.69c  :z=+0.47-1mm 

/  - x=0.69c  :  z=+0.47+1mm 

[a]  Axial  Ax  =  ±lmm. 


[b]  Transverse  Az  =  ±lmm. 


Figure  C.4  Effect  of  a  variation  in  probe  position  on  <j>\  convex  side  (M=2.8) 


Appendix  D.  Mach  5  Wind  Tunnel 


D.l  Overview 

The  AFIT  Mach  5  wind  tunnel  is  a  blow-down  design  (as  illustrated  in  Figure  D.l) 
which  operates  over  a  range  of  Reynolds  numbers.  Since  the  static  temperature  for  Mach  5 
flow  at  standard  day  total  temperature  is  below  the  liquefaction  temperature  of  oxygen, 
the  air  is  heated  by  a  refractive  pebble-bed  heating  system.  The  nominal  chamber  total 
pressure  and  temperature  are  20-32atm  and  350-375K  respectively,  yielding  a  freestream 
Reynolds  number  Re/^«32-75xl06m-1.  The  nominal  tunnel  conditions  experienced  in 
the  tunnel  for  the  present  research  are  shown  in  Table  D.l.  Although  the  tunnel  is  capable 
of  sustaining  longer  run  times,  it  is  typically  run  for  only  10-15  seconds.  The  design, 
construction,  and  instrumentation  of  this  facility  comprised  a  significant  portion  of  this 
research  effort. 

D.2  Air  Supply 

High  pressure  air  is  stored  in  a  175atm  tank  (Figure  D.2)  having  a  capacity  of  1.25  cu¬ 
bic  meters,  and  is  controlled  with  a  high  volume  dome  regulated  control  valve  (Figure  D.3). 
The  tank  is  recharged  with  a  four- stage  compressor  and  dryer  system  (Figure  D.4). 

This  system  provides  enough  air  to  run  the  tunnel  at  a  plenum  pressure  of  20atm  for  over 
2  minutes,  and  requires  about  5  hours  to  re-charge,  although  actual  run-times  are  never 
expected  to  exceed  30  seconds  in  order  to  allow  several  runs  during  a  day. 

Table  D.l  Typical  Mach  5  tunnel  freestream  conditions  (used  for  viscous  CFD). 


Condition 

Value 

M 

4.87 

Pt 

3.23xl06Pa  =  32.0atm 

P 

7.017kPa  =  0.070atm 

Tt 

354K 

T 

61.5K 

u 

766™ 

p 

0.3977^- 

D-l 


Nozzle  Test 
Section 


Diffuser 


Silencer 


S' 


Elec,  heaters  (2)  Regulator 

ss=>-Q--D- 


House  Air 


On/Off 

valve 


Regulator 

gti— 

I 


2500  psi  air 


Figure  D.l  Mach  5  wind  tunnel  schematic 


©GOG 

Dryer 

Compressor 

In  the  present  configuration,  both  the  run  time  and  down  time  are  determined  by 
the  heating  system.  The  system  passes  the  high-pressure  air  through  a  bank  of  refractory 
pebble-bed  heaters  (Figure  D.5)  before  it  enters  the  tunnel.  To  accomplish  the  heater 
design,  computer  programs  (in  both  FORTRAN  and  Mathcad)  were  developed  which  im¬ 
plement  the  heater  design  method  of  Pope  and  Goin.  t86)  Each  of  the  pebble  beds  are 
designed  to  heat  the  air  to  375K  before  entering  the  tunnel,  which  keeps  the  static  tem¬ 
perature  in  the  tunnel  well  above  the  liquefaction  temperature.  This  temperature  can  be 
maintained  or  exceeded  for  several  seconds,  and  is  monitored  by  a  thermocouple  in  the 
settling  chamber.  The  air  can  easily  be  heated  to  higher  or  lower  temperatures  by  heating 
the  pebble  beds  to  other  initial  temperatures.  The  heater  system  is  designed  to  be  very 
flexible,  allowing  the  use  of  single  or  multiple  heaters  during  a  run.  This  is  accomplished  by 
channeling  the  air  from  the  heaters  through  a  manifold  before  entering  the  tunnel,  which 
allows  them  to  be  run  independently  or  in  parallel.  In  practice,  two  heaters  are  usually  run 
in  parallel,  providing  enough  heat  for  three  runs,  each  sustaining  acceptable  temperatures 
for  8-12  seconds  of  10-15  second  runs. 
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Style: 

ffigh  Volume  /  Dome  Regulated  204atm 

Regulator: 

Grove  Powreactor  Dome  Regulator  #200/300 

Manufacturer: 

Eagle  Compressors  Inc., 

Pleasant  Garden,  NC 

Model: 

A.C.P.  Custom  Air  Control  Panel  #G07C0177 

Figure  D.3  Air  control  panel. 
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Style: 

air  cooled  four  stage  compressor 

Power: 

30kW 

Separator: 

mechanical 

Dryer: 

regenerative 

Delivery: 

lm3/min  @  41atm  (free  air  delivery) 

Manufacturer: 

Eagle  Compressors  Inc., 

Pleasant  Garden,  NC 

Model: 

HORIZONTAL  4S50  MKII  #HW40HH3 

Figure  D.4  Compressor  and  dryer. 
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Figure  D.5  Mach.  5  wind  tunnel  pebble  bed  heating  system. 


The  heaters  are  re-heated  by  convection  with  two  heavy-duty  electric  air  heaters  (Re¬ 
heat  Co.  Inc.  Model  HDA-2-12-2).  These  heaters  each  produce  a  minimum  l.okW  when 
powered  by  standard  240  Volt  outlets.  With  both  electric  heaters  in  operation,  it  takes 
approximately  1  hour  to  re-heat  each  pebble-bed  heater  from  an  end-of-run  temperature 
back  to  operating  temperature. 

D.3  Nozzle  &  Test  Section 

The  nozzle  profile  was  designed  using  the  method  of  characteristics  with  modifications 
for  viscous  effects.  The  test  section  of  the  wind  tunnel  has  cross-sectional  dimensions  of 
7.62cm  x  7.62cm  and  a  length  of  26.67cm.  The  test  section  and  nozzle  (Figures  D.6 
and  D.7)  are  milled  from  347  stainless  steel.  The  test  section  side  walls  are  fitted  with 
optical  grade  glass  windows  to  alow  for  photographic  visualization  of  the  flow  (Figure 
D.$).  The  tunnel  ceiling  (Figure  D.9)  is  designed  to  accept  the  WAF  model,  and  the  floor 
(Figure  D.10)  is  fitted  with  slots  for  probe  insertion. 
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Figure  D.7  Wind  tunnel  nozzle  and  test  section.  Shown  with  nozzle  wall  removed  for 
visibility  and  window  replaced  by  aluminum  plug  for  initial  tunnel  calibration. 


D-4  Diffuser  &  Silencer 

An  adjustable  throat-area  diffuser  was  employed,  allowing  it  to  be  tuned  for  optimal 
performance  after  the  tunnel  was  assembled.  A  sound  muffler  (seen  in  Figure  D.2)  which 
is  attached  directly  to  the  diffuser  exit  directs  the  air  back  outside  the  building. 

D.5  Plenum  Chamber  Instrumentation 

Instrumentation  for  the  tunnel  conditions  include  an  upstream  pitot  tube  and  ther¬ 
mocouple  well  ahead  of  the  throat  to  measure  Pi ^  and  Tt0 0.  A  34atm  Endevco  pressure 
transducer  was  used  to  sense  plenum  pressure,  which  is  displayed  and  digitized  by  an  En¬ 
devco  (model  4428A)  signal  conditioner.  Plenum  temperature  is  sensed  with  an  Omega 
K-type  (model  KAIN-18-U-12)  lance  thermocouple,  then  displayed  on  an  Omega  model 
DP41-TC-A  temperature  display.  This  temperature  data  was  used  in  the  data  reduction 
as  well  as  to  ensure  that  the  air  had  been  adequately  heated.  Pressure  and  temperature  is 
sent  to  the  data  acquisition  unit  (Section  3.6). 
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Figure  D.8  Wind  tunnel  test  section  wall  with  windows  (dimensions  in  inches). 

D.6  TSI 3-DOF  Traverse  System 

A  TSI  3-DOF  traverse  system  was  used  for  the  Mach  5  experiments.  Although  this 
traverse  was  equipped  with  its  own  manual  control  center,  safety  concerns  dictated  that 
the  traverse  be  controlled  from  a  remote  location  during  tunnel  operation.  The  control 
software  also  provides  the  ability  to  move  specific  distances  at  prescribed  speeds. 
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0.6630 


Figure  D.9  Wind  tunnel  ceiling  with  model  insert  (dimensions  in  inches). 


Figure  D.ll  Measured  pitot  pressure  at  upstream  location. 

D.  7  Tunnel  Calibration 

When  operated  near  the  conditions  given  in  Table  D.l,  the  tunnel  flow  is  very  two- 
dimensional  over  a  significant  portion  of  the  tunnel  width.  This  is  demonstrated  in  Fig¬ 
ure  D.ll  which  shows  pitot  pressures  measured  at  an  upstream  (£=-9.64r=-7.53c)  loca¬ 
tion  over  a  range  of  span-wise  locations. 


D.8  Numerical  Validation 

The  AFIT  Mach  5  wind  tunnel  was  simulated  numerically  to  provide  some  reassur¬ 
ance  of  the  tunnel  design,  as  well  as  to  gain  some  early  experience  with  the  methods  used 
for  the  wrap-around  fin  simulations.  The  simulations  were  conducted  using  the  General 
Aerodynamic  Simulation  Program  (GASP)  Version  2.CK72h  The  computational  require¬ 
ments  of  these  simulations  were  also  used  to  estimate  the  computational  requirements  of 
the  WAF  simulations  presented  in  this  document.  Several  2-D  and  3-D  solutions  were 
obtained,  the  details  of  which  were  reported  in  the  Prospectus  for  this  research.  *1073 
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Appendix  E.  GASP  Input  Deck  for  M=2.8  Viscous  Simulation 

The  GASP  input  deck  (version  3.0)  for  the  viscous  simulation  at  Mach  2.8  discussed  in 
Chapter  4  is  presented  below.  As  discussed  in  Chapter  4,  a  2-Dimensional  inflow  boundary 
condition  (shown  in  Figure  E.l)  is  specified  on  the  upstream  edge  of  the  first  zone.  Parallel 
flow  and  constant  pressure  through  the  boundary  layer  was  assumed.  This  two-dimensional 
PNS  region  allowed  the  boundary  layer  to  develop  into  a  fully  turbulent  profile  (also  shown 
in  Figure  E.l)  upstream  of  the  blended  body  region. 


Figure  E.l  Space  marched  zone  inflow  boundary  condition  for  and  exit. 
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ieq 


E-ll 


icond 

2 

Vel/Mach  cx 

cy 

cz 

temp/press 

turbi 

tkelref 

2.9  1 

0 

0 

109.6 

0.01 

0.001 

rho_spec 

0.2039 

PHYSICAL  MODEL  #2: 

RANS, 

Viscous 

-  Laminar, 

3rd  order 

CHEMISTRY  &  THERMODYNAMICS 


nspec  nnev  prefDiss  vibRelax  itberm  chemmod 

1  0  0  0  4  Perfect  Gas' 


INVISCID  FLUXES 


imarch 

0 

invflxi  invflxj  invflxk 

2  2  2 
rkapi  rkapj  rkapk  sdm2  sdm4 

0.3333  0.3333  0.3333  0.0000  0.0000 

limi  limj  limk  rk_ven 

222  1.000 


VISCOUS  FLUXES 


isViscous 


1 

visflxi  visflxj  visflxk 


-1 

- 

■1 

-1 

modlmu 

modlk 

imodld 

ivac 

2 

2 

1 

2 

prl 

prt 

scl 

set 

0.72 

0.90 

1.00  0 

.50 

ikeps 

kemin 

wallfunc 

igb 

0 

0 

0 

0 

INITIAL  CONDITIONS 


icond 

2 


Vel/Mach 

CX 

cy 

cz 

temp/press 

turbi  tkelref 

2.9 

1 

0 

0 

109.6 

0.01  0.001 

rho_spec 

0.2039 


REFERENCE  STATE  FOR  B.C.’s 


icond 


ieq 


E-12 


2 

Vel/Mach 

cx 

cy 

cz 

temp/press 

turbi 

tkelref 

2.9 

1 

0 

0 

109.6 

0.01 

0.001 

rho_spec 

0.2039 

******** t****************************************************** 


PHYSICAL  MODEL  #3:  RANS,  Viscous  -  B-L,  3rd  order 


CHEMISTRY  &  THERMODYNAMICS 


uspec  nnev  prefDiss  vibRelax  itherm  chemmod 

100  0  4  ’Perfect  Gas’  1 


INVISCID  FLUXES 


imarch 

0 


invflxi 

invflxj 

invflxk 

2 

2 

2 

rkapi 

rkapj 

rkapk  sdm2 

sdm4 

0.3333 

0.3333 

0.3333  0.0000 

0.0000 

limi  limj  limk 

rk_ven 

2 

2  2 

1.000 

VISCOUS  FLUXES 


isViscous 

2 


visflxi  visflxj  visflxk 


-1 

- 

•1 

-1 

modlmu 

modlk 

imodld 

ivac 

2 

2 

1 

2 

prl 

prt 

scl 

set 

0.72 

0.90 

1.00  0 

'.50 

ikeps 

kemin 

wallfunc 

igb 

0 

0 

0 

0 

INITIAL  CONDITIONS 


icond 

2 


Vel/Mach 

2.9 

rho_spec 

0.2039 

cx 

1 

cy 

0 

cz 

0 

temp/press 

109.6 

turbi 

0.01 

tkelref 

0.001 

REFERENCE 

STATE  FOR  B.C. ’s 

icond 

2 

Vel/Mach 

cx 

cy 

cz 

temp/press 

turbi 

tkelref 

2.9  10  0  109.6  0.01  0.001 

rho.spec 
0.2039 


BLOCK  INFO 


BLOCK  #1:  march  through  zone  1  —  (starting  at  experimental  zb)  (1-60) 
imcont  nswp  ncycle  nwres  mstage  rtolr  rtola 

1  1  250  50  1  1.00e-03  1.00e-10 

mgstyle  nitfine  nitcrct  nitsmth  smthfac  dtfac 

0110  3.00e-01  1.00e+00 

SWEEP  #1:  march  through  zone  3 

nz  iseq  npm  partStyle  iswpdir  iplstt  iplend 
13  1  1  1  1  60 

impl  itmstep  keslv  ichemslv  nplane  inner  mxin  tolin 
2  1  1  1  1  0  5  0.01 

dtmin  dtmax  irelu  nremax  tolreu 

-1.00e+00  -1.00e+00  1  10  1.00e-01 


BLOCK  #2:  march  through  zone  2  (1-100) 

imcont  nswp  ncycle  nwres  mstage  rtolr  rtola 

1  1  250  50  1  1.00e-03  1.00e-10 

mgstyle  nitfine  nitcrct  nitsmth  smthfac  dtfac 

0110  3.00e-01  1.00e+00 

SWEEP  #1:  march  through  zone  4 

nz  iseq  npm  partStyle  iswpdir  iplstt  iplend 
2  3  1  1  1  1  100 

impl  itmstep  keslv  ichemslv  nplane  inner  mxin  tolin 
2  1  1  1  1  0  5  0.01 

dtmin  dtmax  irelu  nremax  tolreu 

-1.00e+00  -1.00e+00  1  10  1.00e-01 


BLOCK  #3:  march  through  zone  3  (1-80) 

imcont  nswp  ncycle  nwres  mstage  rtolr  rtola 

1  1  250  50  1  1.00e-03  1.00e-10 

mgstyle  nitfine  nitcrct  nitsmth  smthfac  dtfac 

0110  3.00e-01  1.00e+00 

SWEEP  #1:  march  through  zone  5 

nz  iseq  npm  partStyle  iswpdir  iplstt  iplend 
3  3  1  1  1  1  80 

impl  itmstep  keslv  ichemslv  nplane  inner  mxin  tolin 
2  1  1  1  1  0  5  0.01 

dtmin  dtmax  irelu  nremax  tolreu 

-1.00e+00  -1.00e+00  1  10  1.00e-01 


BLOCK  #4:  march  through  zone  4 

imcont  nswp  ncycle  nwres  mstage  rtolr  rtola 

1  1  250  50  1  1.00e-03  l.OOe-10 

mgstyle  nitfine  nitcrct  nitsmth  smthfac  dtfac 

0110  3 . 00e-01  1.00e+00 

SWEEP  #1:  march  through  zone  6 


E-14 


nz  iseq  npm  part Style  iswpdir  iplstt  iplend 
4  3  1  1  1  1  24 

impl  itmstep  keslv  ichemslv  nplane  inner  mxin  tolin 
2  1  1  1  1  0  5  0.01 

dtmin  dtmax  irelu  nremax  tolreu 

-5.00e-01  -5.00e-01  1  10  l.OOe-Ol 


BLOCK  #5:  global  by  fin  -  Laminar  -  sequence  3 

imcont  nswp  ncycle  nwres  mstage  rtolr  rtola 

0  8  2000  100  1  1.00e-04  1.00e-10 

mgstyle  nitfine  nitcrct  nitsmth  smthfac  dtfac 

0110  3.00e-01  1.00e+00 

SWEEP  #1: 

nz  iseq  npm  partStyle  iswpdir  iplstt  iplend 

10  3  2  1  2  1  96 

impl  itmstep  keslv  ichemslv'  nplane  inner  mxin  tolin 

4  11  1  115  0.01 

dtmin  dtmax  irelu  nremax  tolreu 

-i.00e-00  -1.00e-00  0  10  1.00e-01 

SWEEP  #2: 

nz  iseq  npm  partStyle  iswpdir  iplstt  iplend 

9  3  2  1  2  1  96 

impl  itmstep  keslv  icbemslv  nplane  inner  mxin  tolin 

4  11  1  115  0.01 

dtmin  dtmax  irelu  nremax  tolreu 

-1.00e-00  -1 .00e-00  0  10  1.00e-01 

SWEEP  #3: 

nz  iseq  npm  partStyle  iswpdir  iplstt  iplend 

5  3  2  1  2  1  40 

impl  itmstep  keslv  ickemslv  nplane  inner  mxin  tolin 

4  1  1  1  1  15  0.01 

dtmin  dtmax  irelu  nremax  tolreu 

-1.00e-00  -1.00e-00  0  10  1.00e-01 

SWEEP  #4: 

nz  iseq  npm  partStyle  iswpdir  iplstt  iplend 

6  3  2  1  2  1  40 

impl  itmstep  keslv  ichemslv  nplane  inner  mxin  tolin 

4  11  1  115  0.01 

dtmin  dtmax  irelu  nremax  tolreu 

-1.00e-00  -1.00e-00  0  10  1.00e-01 

SWEEP  #5: 

nz  iseq  npm  partStyle  iswpdir  iplstt  iplend 

7  3  2  1  2  1  40 

impl  itmstep  keslv  ichemslv  nplane  inner  mxin  tolin 

4  11  1  115  0.01 

dtmin  dtmax  irelu  nremax  tolreu 

-i.00e-00  -1.00e-00  0  10  1.00e-01 

SWEEP  #6: 

nz  iseq  npm  partStyle  iswpdir  iplstt  iplend 

8  3  2  1  2  1  40 

impl  itmstep  keslv  icbemslv  nplane  inner  mxin  tolin 


E-15 


4  11  1  115  0.01 

dtmin  dtmax  irelu  nremax  tolreu 

-l.OOe-OO  -1.00e-00  0  10  1.00e-01 

SWEEP  #7: 

nz  iseq  npm  part Style  iswpdir  iplstt  iplend 

11  3  2  1  3  1  20 

impl  itmstep  keslv  ichemslv  nplane  inner  mxin  tolin 

4  11  1  115  0.01 

dtmin  dtmax  irelu  nremax  tolreu 

-1.00e-00  -1.00e-00  0  10  1.00e-01 

SWEEP  #8: 

nz  iseq  npm  partStyle  iswpdir  iplstt  iplend 

12  3  2  1  3  1  20 

impl  itmstep  keslv  ichemslv  nplane  inner  mxin  tolin 

4  11  1  115  0.01 

dtmin  dtmax  irelu  nremax  tolreu 

-1.00e-00  -1.00e-00  0  10  1.00e-01 


BLOCK  #6:  global  by  fin  -  Laminar  -  sequence  2 

imcont  nswp  ncycle  nwres  mstage  rtolr  rtola 

0  8  3000  100  1  1 . 00e-04  1.00e-10 

mgstyle  nitfine  nitcrct  nitsmth  smthfac  dtfac 

0110  3 .00e-01  1.00e+00 

SWEEP  #1: 

nz  iseq  npm  partStyle  iswpdir  iplstt  iplend 

10  2  2  1  2  1  96 

impl  itmstep  keslv  ichemslv  nplane  inner  mxin  tolin 

4  11  1  115  0.01 

dtmin  dtmax  irelu  nremax  tolreu 

-1.00e-00  -1.00e-00  0  10  1.00e-01 

SWEEP  #2: 

nz  iseq  npm  partStyle  iswpdir  iplstt  iplend 

9  2  2  1  2  1  96 

impl  itmstep  keslv  ichemslv  nplane  inner  mxin  tolin 

4  11  1  115  0.01 

dtmin  dtmax  irelu  nremax  tolreu 

-1.00e-00  -1.00e-00  0  10  1.00e-01 

SWEEP  #3: 

nz  iseq  npm  partStyle  iswpdir  iplstt  iplend 

5  2  2  1  2  1  40 

impl  itmstep  keslv  ichemslv  nplane  inner  mxin  tolin 

4  11  1  115  0.01 

dtmin  dtmax  irelu  nremax  tolreu 

-1.00e-00  -l.OOe-OO  0  10  1.00e-01 

SWEEP  #4: 

nz  iseq  npm  partStyle  iswpdir  iplstt  iplend 

6  2  2  1  2  1  40 

impl  itmstep  keslv  ichemslv  nplane  inner  mxin  tolin 

4  11  1  115  0.01 

dtmin  dtmax  irelu  nremax  tolreu 

-1.00e-00  -1.00e-00  0  10  1.00e-01 


E-16 


SWEEP  #5: 
nz  iseq 

npm 

partStyle 

iswpdir 

iplstt  iplend 

7  2 

2 

1 

2 

1  40 

impl  itmstep 

keslv  ichemslv  nplane  inner  mxin 

4 

1 

1 

1 

115 

dtmin 

dtmax  irelu 

nremax 

tolreu 

-i.00e-00 

-l.OOe- 

-00  0 

10 

1.00e-01 

SWEEP  #6: 
nz  iseq 

npm 

partStyle 

iswpdir 

iplstt  iplend 

8  2 

2 

1 

2 

1  40 

impl  itmstep  keslv  ichemslv  nplane  inner  mxin 


4 

1 

1 

1 

1  1  5 

dtmin 

dtmax 

irelu 

nremax 

tolreu 

-1.00e-00 

-l.OOe- 

■00  0 

10 

i.00e-01 

SWEEP  #7: 

nz  iseq 

npm 

partStyle 

iswpdir 

iplstt  iplend 

11  2 

2 

1 

3 

1  20 

impl  itmstep 

keslv  ichemslv  nplane  inner  mxin 

4 

1 

1 

1 

115 

dtmin 

dtmax 

irelu 

nremax 

tolreu 

-i.00e-00 

-l.OOe- 

■00  0 

10 

1.00e-01 

SWEEP  #8: 

nz  iseq 

npm 

partStyle 

iswpdir 

iplstt  iplend 

12  2 

2 

1 

3 

1  20 

impl  itmstep 

keslv  ichemslv  nplane  inner  mxin 

4 

1 

1 

1 

115 

dtmin 

dtmax 

irelu 

nremax 

tolreu 

-1.00e-00 

-l.OOe- 

■00  0 

10 

1.00e-01 

BLOCK  #7:  global  by  fin  -  Laminar  -  sequence  1 

imcont  nswp  ncycle  nwres  mstage  rtolr  rtola 


0  8  4000  100 

1 

1.00e-03 

1.00e-10 

mgstyle  nitfine 

nitcrct  nitsmth 

smthfac 

dtfac 

0  1 

SWEEP  #1: 

1 

0 

3.00e-01 

1.00e+00 

nz  iseq  npm 

partStyle 

iswpdir 

iplstt 

iplend 

10  1  2 

1 

2 

1 

96 

impl  itmstep 

keslv  ichemslv  nplane  inner  mxin 

4  1 

1 

1 

1  1 

5 

dtmin  dtmax  irelu 

nremax 

tolreu 

-1.00e-00  -l.OOe- 

-00  0 

10 

l.OOe-Ol 

SWEEP  #2: 
nz  iseq  npm 

partStyle 

iswpdir 

iplstt 

iplend 

9  12 

1 

2 

1 

96 

impl  itmstep 

keslv  ichemslv  nplane  inner  mxin 

4  1 

1 

1 

1  1 

5 

dtmin  dtmax  irelu 

nremax 

tolreu 

-1.00e-00  -l.OOe- 

-00  0 

10 

l.OOe-Ol 

SWEEP  #3: 
nz  iseq  npm 

partStyle 

iswpdir 

iplstt 

iplend 

5  12 

1 

2 

1 

40 

E-17 

tolin 

0.01 


tolin 

0.01 


tolin 

0.01 


tolin 

0.01 


tolin 

0.01 


tolin 

0.01 


impl  itmstep  keslv  ichemslv  nplane  inner  mxin  tolin 

4  1  1  1  1  15  0.01 

dtmin  dtmax  irelu  nremax  tolreu 

-1 .OOe-OO  -1.00e-00  0  10  1.00e-01 

SWEEP  #4: 

nz  iseq  npm  partStyle  iswpdir  iplstt  iplend 

6  12  1  2  1  40 

impl  itmstep  keslv  ichemslv  nplane  inner  mxin  tolin 

4  11  1  115  0.01 

dtmin  dtmax  irelu  nremax  tolreu 

-1 . 00e-00  -1.00e-00  0  10  1.00e-01 

SWEEP  #5: 

nz  iseq  npm  partStyle  iswpdir  iplstt  iplend 

7  12  1  2  1  40 

impl  itmstep  keslv  ichemslv  nplane  inner  mxin  tolin 

4  11  1  115  0.01 

dtmin  dtmax  irelu  nremax  tolreu 

-1.00e-00  -1.00e-00  0  10  1.00e-01 

SWEEP  #6: 

nz  iseq  npm  partStyle  iswpdir  iplstt  iplend 

8  12  1  2  1  40 

impl  itmstep  keslv  ichemslv  nplane  inner  mxin  tolin 

4  11  1  115  0.01 

dtmin  dtmax  irelu  nremax  tolreu 

-1.00e-00  “1.00e-00  0  10  1.00e-01 

SWEEP  #7: 

nz  iseq  npm  partStyle  iswpdir  iplstt  iplend 

11  1  2  1  3  1  20 

impl  itmstep  keslv  ichemslv  nplane  inner  mxin  tolin 

4  1  1  1  1  1  50. 01 

dtmin  dtmax  irelu  nremax  tolreu 

-1.00e-00  -1.00e-00  0  10  1.00e-01 

SWEEP  #8: 

nz  iseq  npm  partStyle  iswpdir  iplstt  iplend 

12  1  2  1  3  1  20 

impl  itmstep  keslv  ichemslv  nplane  inner  mxin  tolin 

4  11  1  115  0.01 

dtmin  dtmax  irelu  nremax  tolreu 

-i.00e-00  -1.00e-00  0  10  1.00e-01 


BLOCK  #7:  global  by  fin  -  Turbulent  B-L  -  sequence  1 

imcont  nswp  ncycle  nwres  mstage  rtolr  rtola 

0  8  5000  100  1  1.00e-03  1.00e-10 

mgstyle  nitfine  nitcrct  nitsmth  smthfac  dtfac 

0110  3 .00e-0i  1.00e+00 

SWEEP  #1: 

nz  iseq  npm  partStyle  iswpdir  iplstt  iplend 
10  1  3  1  2  1  96 

impl  itmstep  keslv  ichemslv  nplane  inner  mxin  tolin 
4  11  1  115  0.01 

dtmin  dtmax  irelu  nremax  tolreu 


E-18 


0 


10 


l.OOe-Ol 


-1.00e-00  -1.00e-00 

SWEEP  #2: 

nz  iseq  npm  partStyle  iswpdir  iplstt  iplend 

9  13  1  2  1  96 

impl  itmstep  keslv  ichemslv  nplane  inner  mxin  tolin 

4  11  1  115  0.01 

dtmin  dtmax  irelu  nremax  tolreu 

-1.00e-00  -1.00e-00  0  10  l.OOe-Ol 

SWEEP  #3: 

nz  iseq  npm  partStyle  iswpdir  iplstt  iplend 

5  12  1  2  1  40 

impl  itmstep  keslv  ichemslv  nplane  inner  mxin  tolin 

4  11  1  115  0.01 

dtmin  dtmax  irelu  nremax  tolreu 

-1.00e-00  -1.00e-00  0  10  l.OOe-Ol 

SWEEP  #4: 

nz  iseq  npm  partStyle  iswpdir  iplstt  iplend 

6  12  1  2  1  40 

impl  itmstep  keslv  ichemslv  nplane  inner  mxin  tolin 

4  11  1  115  0.01 

dtmin  dtmax  irelu  nremax  tolreu 

-1.00e-00  -1.00e-00  0  10  l.OOe-Ol 

SWEEP  #5: 

nz  iseq  npm  partStyle  iswpdir  iplstt  iplend 

7  12  1  2  1  40 

impl  itmstep  keslv  ichemslv  nplane  inner  mxin  tolin 

4  11  1  115  0.01 

dtmin  dtmax  irelu  nremax  tolreu 

-1.00e-00  -I . 00e-00  0  10  1.00e-01 

SWEEP  #6: 

nz  iseq  npm  partStyle  iswpdir  iplstt  iplend 

8  12  1  2  1  40 

impl  itmstep  keslv  ichemslv  nplane  inner  mxin  tolin 

4  11  1  115  0.01 

dtmin  dtmax  irelu  nremax  tolreu 

-1.00e-00  -1.00e-00  0  10  l.OOe-Ol 

SWEEP  #7: 

nz  iseq  npm  partStyle  iswpdir  iplstt  iplend 

11  1  3  1  3  1  20 

impl  itmstep  keslv  ichemslv  nplane  inner  mxin  tolin 

4  11  1  115  0.01 

dtmin  dtmax  irelu  nremax  tolreu 

-1.00e-00  -1.00e-00  0  10  l.OOe-Ol 

SWEEP  #8: 

nz  iseq  npm  partStyle  iswpdir  iplstt  iplend 

12  1  3  1  3  1  20 

impl  itmstep  keslv  ichemslv  nplane  inner  mxin  tolin 

4  11  1  115  0.01 

dtmin  dtmax  irelu  nremax  tolreu 

-1.00e-00  -1.00e-00  0  10  l.OOe-Ol 


E-19 


